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Tidal Mixing in the Narrows of the Mersey Estuary 


P. Hughes 


(Received 1958 June 2) 


Summary 


Observations of salinity in the Narrows of the River Mersey have 
shown that vertical mixing due to tidally induced turbulence may be 
complete at mid-tide. As tidal currents abate horizontal differential 
velocities, which constitute a seaward flow in the upper layer and an 
upriver flow near the bottom, become apparent and produce vertical 
salinity gradients. A second series of observations showed that vertical 
homogeneity may not be reached at any stage of the tide if the riverflow 
is high. Applying the theory for a vertically mixed estuary, coefficients 
of longitudinal eddy diffusion were computed for a low riverflow of 
26 m3/s and gave a mean value of 160m?/s. During a period of greater 
fresh water influx (103 m3/s) when conditions were less clearly defined, 
the average coefficient was considerably higher at 360 m2/s. 


1. Introduction 


The part of the Estuary of the River Mersey called the Narrows is a straight 
deep channel, with a fairly regular bottom, which connects Liverpool Bay with 
the upper basin. From New Brighton the Narrows extends 9-7km (6ml) up- 
stream and is only 1-2km wide at its narrowest, although in places at low water 
the depth exceeds 21-3 m (70ft). The estuary then widens into the upper basin 
extending as far as Howley weir, Warrington, above which the river is non-tidal 
except for very high tides. The upper basin has a maximum width of about 4°8 
km (3 ml) but it is shallow, and at low water many sandbanks can be seen and the 
water in the channels is usually only a few metres deep. 

For a tide reaching 9-4m (31 ft) above the Liverpool Bay Datum at Princes 
Stage, Liverpool, the intertidal capacity of the estuary above Rock Light, New 
Brighton, is about 550 x 10m? whereas at low water the capacity is only about 
150 x10®m%. The high-water capacity from the Rock Light, New Brighton, to a 
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section 8km upriver is 195-4 x 106m, whereas the low-water capacity is 99 x 106 m3. 
Therefore the mean volume of water in this region is-147-2 x 10 m3 and the inter- 
tidal volume is 96.4 x 106 m3. From these data the mean current over a 9-4 m flood 
or ebb tide may be computed. Considering first the largest cross-section at Rock 
Light, which has a mean area of 24-7 10°m2, the mean current is 1-om/s, and assum- 
ing a sinusoidal variation of tidal current the maximum current is 1-6 m/s. Simi- 
larly for the smallest cross-section, which is situated 4km from Rock Light and 
has a mean area of 15 x 10% m2, the mean current is estimated at 1-4 m/s and the 
maximum current at 2:2 m/s. 

The riverflow due to the influx of fresh water from the run-off areas is small 
compared with the large tidal currents. For the first series of observations in 
1956 May the flow amounted to 25-7 m3/s whereas between 1957 March 17 and 23 
it varied from 205 to 68 m/s. 

The observational programme described here was by way of a reconnaissance, 
with a view to a more detailed study of mixing processes. The Narrows was 
chosen as the shape of the channel is simple, whereas the rest of the estuary presents 
a complicated system of varying width and depth. Moreover, being a small 
area it was possible to obtain a detailed picture of the variables concerned in a short 
time. 


2. Observations of temperature and salinity 


Sampling of the estuarine water was carried out from a survey launch provided 
by the Mersey Docks and Harbour Board, between 1956 May 22 and 1956 May 28, 
and again on 1957 March 21 and 22. The approximate positions of sampling are 
shown in Figure 1. At each station samples were obtained from the surface, mid- 
depth and 1:2m (4ft) above the bottom. The surface sample was taken with a 
canvas bucket, while the other two were collected with a Nansen—Pettersson 
insulating water bottle. In the first series, one cross-section of the river was covered 
each day by sampling at three stations, one fairly close to each bank and the other 
in the middle. A station was occupied every half hour, so that surface, interme- 
diate and bottom samples were obtained every hour in the middle and every two 
hours at each of the outer stations. The operation was repeated on different days. 
For sections I and II, ebb tides only were observed, at III a whole tidal cycle 
was sampled, while at IV and V observations were carried out over flood tides. 
The period was chosen because the predicted variation of tidal range was small. 
The changes from one tidal cycle to the next were observed by a series of surface 
samples taken at New Brighton pier. 

The results indicated that for the second series, taken during a period of fairly 
high fresh-water influx in 1957 March, it was sufficient to sample at a mid-river 
station at each end of the Narrows. 

Temperatures were measured with mercury-in-glass thermometers and are 
correct to o-1 °C. The salinity of each sample was measured in the laboratory by 
an electrical conductivity method. The apparatus consisted of a high quality 
resistance bridge with a 600 c/s tone source and an earphone detector. The sample 
was contained in a glass cell and the resistance between two platinum electrodes 
in this cell constituted the unknown arm of the bridge. Calibration was effected 
by using diluted sea-water of known salinity. The electrical conductivity is sensi- 
tive to temperature changes and it was necessary to maintain the temperature of 
the cell to within 0-o1°C in order to give an accuracy of 0-05 per cent in the 
salinity. 
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Fic. 2.—Variation of salinity with time at mid-river stations 


(b) Station 5, 1956 May 23, 
(d) Station 11, 1956 May 25. 
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Fig. 2. Contd.—Variation of salinity with time at mid-river stations 
(e) Station 16, 1957 March 21, (f) Station 17, 1957 March 22. 
— @© — surface, —O— intermediate, — A — bottom. 


3- The distribution of salinity, temperature and density 


The salinity distribution during the period of the first series of observations in 
1956 May will be considered first. The variation of salinity with time and depth 
at each mid-river station is shown in Figure 2(a) to (d). The increase in salinity 
at all levels as the tide floods and its decrease as the tide ebbs is obvious. The 
difference between the extreme values over a tidal cycle is about 4 parts per thou- 
sand (°/,.). The kink in the curve of surface salinity at Station 8 on May 24 
was probably due to the local influx of fresh water across the surface during a rain 
shower. The variation of salinity with depth was small or absent when the tidal 
current was strong, so that for the mid-tide period the isohalines were sensibly 
vertical. Near high and low water, however, there were significant salinity grad- 
ients, with the salinity increasing from the surface to the bottom. 

Because of the small number of observations, only a qualitative study of the 
difference of surface salinity across the Narrows was possible. It shows on the 
average a small increase from the Liverpool to the Birkenhead side of the river. 
The concentration of fresher, lighter water on the right-hand side of a wide estuary, 
looking in the direction of the residual flow, has been observed previously and is 
probably due to the Earth’s rotation (Pritchard 1955). In the narrow channel 
under consideration, however, this phenomenon may be due to local topo- 
graphical effects and there are too few observations on which to base a definite 
conclusion. 

An attempt has been made to present a synoptic picture of the salinity distribu- 
tion in a vertical section along the middle of the Narrows at the time of high and 
low water (registered on the Princes Stage tide gauge). Corrections were applied 
for the change in conditions over the period of observations by referring to the 
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surface salinity at New Brighton. Despite a few irregularities, the isohaline 
diagrams (Figures 3 and 4) give a clear indication of stratification. 
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Fic. 3.—Isohalines at high water in a vertical section along the middle of 
the Narrows (distance in km from Rock Light, New Brighton). 


The salinity differences between surface and bottom were not at a maximum at 
high or low water, but continued to increase to about 1 °/,, some 40 min after 
high water, and to more than 1 °/,, over one hour after low water. The time 
differences given here are averages of the available observations, and typical 


variations of salinity with time between surface and bottom are given in 
Figure 5. 


-6 -4 -3 -2 -1 © 
Distance (km). 


Fic. 4.—Isohalines at low water in a vertical section along the middle 
of the Narrows (distance in km from Rock Light, New Brighton). 


From the salinity data for each cross-section, a mean value over a tidal cycle 
was obtained. Where observations extended over part of a tidal cycle extrapolation 
was necessary. The observations at New Brighton were used again to correct 
for the variation in the tides and meteorological conditions during the sampling 
period. The final mean salinities S are shown in Figure 6, as a function of distance 
along the Narrows. The main point of interest is that, excluding the innermost 
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Fic. 5.—Surface (S,) and bottom salinity (S;) differences. 
(S, > S, then difference > 0). 


(a) 1956 May 22, H.W. o922, (b) 1956 May 25, L.W. 0553. 


section, the salinity gradient along the limited portion of the estuary was sensibly 
constant. The salinity distribution is given by 


S = 30°11 +0°236x 
where x is measured in kilometres from Rock Light, New Brighton, positive to 
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Fic. 6. Distribution of mean salinity along the Narrows (distance in km 
from Rock Light, New Brighton). 


(a) 1956 May, (b) 1957 March. 


seaward. The mean fractional fresh-water concentration, f, over a section, may 
be determined from the following relation: 
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The source sea-water, So, was taken as the mean surface salinity during the period 
of observations at the Liverpool Bar Light Vessel, which is stationed about 22:5 
km (14 naut. ml) outside the estuary. So was found to be 31-68 °/,, and thus 


= 0°049—0°007449x. 


The observations made on 1957 March 21 and 22, when the riverflow was 
substantially higher, will now be considered. Only two sampling stations were 
used and a linear distribution of salinity along the Narrows was assumed. As 
before, corrections for non-stationary conditions were applied to the mean salinity. 
The results are shown in Figure 2 (e) and (f), and indicate a steeper mean horizontal 
salinity gradient and lower overall values of the mean salinity. At each of the 
two mid-river stations a vertical salinity gradient was always apparent and there 
was thus no period when the isohalines were vertical. The differences in salinity 
between surface and bottom at slack water were larger than obtained in the pre- 
vious series of observations. The salinity distribution is given by 


S = 27°364+0°375x. 


The variations of temperature were comparatively small, but where a salinity 
gradient existed a decrease in the temperature from the surface to the bottom was 
generally observed. At any station the temperature decreased on the flood tide 
and increased again on the ebb tide, with a total range of little more than 1°C. 
The vertical temperature gradient at high water was of the order of 0-02° C/m. 


Following the normal oceanographic practice the density is expressed in the 
contracted form, o;, defined by the equation 


ot 
a, (: + g/cm? 
1 000 
where p,; is the density at atmospheric pressure for a salinity S°/,, and a tempera- 
ture ¢°C. This dependence of o; on the salinity and the temperature has been 
determined experimentally and is tabulated (Knudsen 1901; Mathews 1932). 
The distribution of o, in a vertical section along the middle of the Narrows at 
high and low water, which has been computed from the temperature and salinity 


data, is shown in Figures 7 and 8. There is a close resemblance between the distri- 
butions of the salinity and the density. 


4- Computation of the apparent longitudinal eddy diffusivity 


Stommel’s theory of mixing in a vertically homogeneous tidal estuary has the 
advantage that no assumption is made about the physical nature of the longitudinal 
diffusivity, which is deduced from the observed distribution of fresh water (Stom- 
mel 1953). To a certain extent the diffusion of a pollutant in a tidal estuary may 
be similar to its diffusion in turbulent pipe flow, for which a theory has recently 
been given by Taylor (1954). He has shown that the apparent longitudinal com- 
ponent of turbulent diffusion depends on a balance between the transverse diffusion 
and the advection due to the variation of the mean velocity across a section. 
The theory has been extended by Ellison (1957) to the case of a broad open channel, 
but the mean velocity profile and the vertical eddy diffusivity must be known in 
order to obtain numerical values of the longitudinal diffusivity. Moreover in these 
treatments the advection is due to the variations in the velocity caused by the 
friction at the boundary, but where there is a density gradient a further source of 
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advection exists which may be more important. In the estuary under consideration 
little is known of these quantities, and it is impossible at the present time therefore 
to make more than a qualitative use of the ideas. The data have been analysed on 
the lines of Stommel’s theory to give empirical values of the apparent longitudinal 
diffusivity, it being understood that these values may include the effects of the 
physical processes discussed. 
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Fic. 7.—The distribution of the density, o;, at high water in a vertical 
section along the middle of the Narrows. (Distance in km from Rock 
Light, New Brighton). 
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Fic. 8.—The distribution of the density, o;, at low water in a vertical 
section along the middle of the Narrows. (Distance in km from Rock 
Light, New Brighton.) 


In Stommel’s theory of tidal mixing it is assumed that (a) changes in salinity 
both vertically and horizontally across the estuary are small compared with those 
along the river and (b) conditions do not change from one tidal cycle to the next. 

The net flux of a pollutant, F(x) depends on the transport downstream by the 
river flow R, and the effective turbulent diffusion upstream due to tidal mixing. 
This is expressed in the equation 


F(x) =R Ka“ 
=, 
(x 


where c is the concentration of pollutant, A is the cross-sectional area of the section 
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under consideration and K is the longitudinal coefficient of eddy diffusion. If 
the fresh-water content is considered as the pollutant then 


F(x) = R = Rf- kal. 


and 


A(dfjdx)’ 
from which K may be computed. 

The estimation of riverflow followed the method described in the Water 
Pollution Research Technical Paper No. 7 (H.M.S.O. 1938), where the flow from 
each of the streams making a contribution to the total influx of fresh water into 
the Estuary was considered. When gaugings were not available, estimates of 
flow based on the size and character of the drainage area were used. During the 
first period of observations the riverflow was small with an average of 25-7 m3/s. 
On the other hand, in 1957 March the riverflow was changing fairly rapidly and 
dropped from 205 m3/s on March 17, to 68 m3/s on March 23. In view of the change 
in conditions K for this period was computed from two values of the riverflow R, 
namely 103-1 m/s, which is the average for the two days of observations, and 151°8 
m/s, the average over the period March 16 to March 22. 

The mean cross-sectional area, over a tide, for any position was calculated from 
the cross-sectional areas at high water and low water on a tide reaching 31 ft at 
Princes Stage. This information was made available by the Marine Surveyor 
(Mersey Docks & Harbour Board). 

The computed values of the eddy diffusivity are given in Table 1. 


Table 1 
Computation of eddy diffusivity 
Distance Cross- Kn (m*/s) 
Date from New _ sectional 
Brighton area 
(km) 10% m? 
(R = 25-69 (m*/s)) 
1956 May 24°7 132°8 
19°4 165-4 
—4°85 17°0 185°5 
—8-0 160°8 
Mean K 161'0 
(i) (ii) 
(R = 103-05 (m*/s)) (R = 151-78 (m*/s)) 
1957 March —2-4 19°7 368-1 542°2 
19°15 349°6 515°0 
Mean K 358°9 528-6 


For similar conditions of low riverflow the computed values of the coefficient 
of eddy diffusion are of the same order as those obtained by Stommedl in his analysis 
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of Bassindale’s River Severn salinity data. Stommel’s values for K varied widely, 
and an average between his sections 13 and 20, where the fresh-water content was 
roughly the same as in the Narrows for the May sampling period, was 63-8 m2/s. 


5. An estimate of the differential flow between surface and bottom 


There is a tendency for fresher water to move seawards in an upper layer with 
an upriver flow of more saline water near the bottom. An estimate of the lower limit 
of this differential velocity may be made with reference to the high-water isohaline 
diagram and the curves showing the variation of salinity at different depths during 
a tidal cycle. On 1956 May 25 the salinity gradient was absent until an hour and 
ahalf before high water. From then until high water the differential flow increased 
the salinity gradient, despite the effect of turbulent mixing, until at high water a 
distance of 2 km separated the extremes of the 30-00 °/,, isohaline. This implies a 
mean differential current between surface and bottom during the period of 37cm/s. 
The action of turbulence before high water would tend to reduce the salinity 
gradient, so that the above estimate represents a lower limit. It must be 


noticed, however, that over a complete tidal cycle the mean current difference 
will be less. 


6. Flushing time 
The average length of time required for river water to move through a portion 
of an estuary is known as the flushing time 7, and is given by the formula: 
T= Q')R, 


where Q’ is the total amount of river water accummulated in the region considered, 
in our case the Narrows, and R is the riverflow. 
Results for the Narrows are set out in Table 2. 


Table 2 
Mean Mean vol. of Q’ Mean vol. of R Y i 
f water in accumulated 
Narrows fresh water 
1956 m? m? m3/d (days) 
May 22-25 00795 147°2 X10° 2°22 X10° 5°3 
1957 
March 21-22 0:2281 147°2 X 10° 33°58 X 10° 8-91 X 10° 38 \ 
13°12 X10° 26 


For a low river flow it is seen that the average length of time for river water to 
move through the region is 5-3 days, which agrees fairly well with a value of 6-2 
days given in the Water Pollution Research Technical Paper No. 7 (1938) for 
similar conditions. As might be expected, a larger river flow reduces the flushing 


time. 
7. Discussion 


Tidal currents in the Narrows of the River Mersey are mainly responsible 
for the turbulent mixing which takes place, while the seaward transport due to 
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the riverflow is comparatively small. One would expect therefore that the coeffi- 
cients of eddy diffusivity computed from data obtained in 1956 May would charac- 
terize the state of mixing for most riverflows. The results of later observations 
do not support this conclusion, however, as the computed coefficients for 1957 
March 20 and 21 are somewhat higher. The answer is perhaps that conditions in 
the Narrows did not strictly satisfy the initial assumptions of Stommel’s theory, 
insisting as it does on stationary conditions and vertical homogeneity. The 
riverflow was lessening during the later set of observations while the samples 
taken at New Brighton suggest that the salinity was still adjusted to the higher 
riverflow. As regards the vertical distribution of salinity, it can be seen that during 
the two day sampling period the water did not become vertically homogeneous, 
even at times of maximum horizontal tidal current. The effects of tidal mixing 
were still clear, however, as the large difference in salinity between surface and 
bottom, which occurred near low water, became smaller as the tidal currents in- 
creased. This decrease continued until about one and a half hours before high 
water in one case, and about two and a half hours in the other, when minimum 
differences were observed. 

In an estuary which has a source of fresh water at one end and pure sea-water 
at the other, the absence of stratification indicates intense vertical turbulence. 
If there is very little turbulence the motion produces a two-layer system. Vertical 
mixing raises the potential energy of a column of stratified water and the energy is 
derived from the kinetic energy of the turbulence or from the mean motion 
(Proudman 1953). Differential movements between surface and bottom always 
tend to degenerate the distribution to the two-layer type with a dense layer below 
a fresh outflow. The observed distribution of salt is indicative of the balance be- 
tween the turbulence and this horizontal motion. 

Strong vertical turbulence accompanies the tidal current in the Narrows and 
may completely destroy any vertical salinity gradients during the part of the tidal 
cycle when horizontal tidal currents are swift. The condition depends to a certain 
extent on the influx of fresh water. High fresh water run-off may result in strati- 
fication being maintained throughout the tidal cycle, as was observed in 1957 
March. It must be remembered that the somewhat lower tidal ranges would also 
have an effect. 

The presence of stratification and the form of the distribution of salinity, with 
the isohalines sloping at high and low water but vertical or near vertical at mid- 
tide, show that there is a horizontal differential flow between the surface and 
bottom layers. The estimate of this current difference in Section 5 indicates that 
it can exceed 37 cm/s even at a time of low fresh-water influx, while the average 
over a whole tidal cycle will be less. 

Pritchard has classified estuarine circulation patterns into four general types 
(A, B, C and D). Type A is the highly stratified, salt wedge circulation, type B 
exhibits moderate stratification, having a deep saline inflow and a fresher outflow 
nearer the surface, while types C and D are vertically homogeneous, but differ 
from each other in that type C shows marked salinity differences across the estuary 
and must therefore be fairly wide, whereas a type D estuary has salinity changes 
only along its length. When considered over one or more tidal cycles the Narrows 
exhibits type B characteristics. During a tidal cycle, however, the turbulent mixing 
in this part of the estuary may be so strong that for a period it becomes of type D. 
A low riverflow aids this transition, whereas a high fresh-water influx and low tidal 
ranges help to maintain it as a type B estuary. 
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An extensive programme of sampling for various tidal ranges and riverflows 
would be necessary to give a full account of the effects of different variables. Obser- 
vations of velocity should reveal the extent of the surface and bottom layers and 
give a direct measure of their relative speeds. 
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On the Chemical Composition of the Outer Core* 


Gordon J. F. MacDonald and Leon Knopoff 


(Received 1958 August 6) 


Summary 


The average composition of the Earth is assumed to be the same 
as that of the chondritic meteorites. This assumption is consistent 
with the present rate of heat production in the Earth. A core made 
up solely of iron-nickel is inconsistent with the chondritic model unless 
the ratio of the abundances of (Fe+Mg)/Si is unity in the mantle. 
The condition (Fe+Mg)/Si = 1 requires that a major part of the 
mantle be made of material having the composition of pyroxene or 
pyroxene-garnet (eclogite). If the mantle is assumed to be made up of a 
400 km shell of the composition of eclogite underlain by a material of 
dunite composition, the chondritic model requires a core of composi- 
tion (Fe, Ni);.¢Si. 

The representative atomic number (Z) has been introduced to 
characterize compounds at moderately high pressures. From asymptotic 
relations in the Thomas—Fermi-Dirac model, a basis is formed for the 
computation of Z. Bullen’s density distribution gives Z for the outer 
core of about 22. The representative atomic number of Fe}.¢Si is 22°8. 
The composition of the core deduced from equation of state arguments 
is thus consistent with the chondritic model of the Earth. 

The electrical conductivity of (Fe, Ni);.6Si should be a third to a 
tenth that of pure iron for conditions obtained in the core. Such a 
conductivity is consistent with the requirements of the dynamo theory 
of the magnetic field. 


1. Introduction 


The outer core of the Earth is generally assumed to be a liquid iron-nickel 
alloy. The evidence advanced in support of this view includes (1) the existence of 
iron meteorites, (2) the fact that the density of iron is in general agreement with 
the density of the core as estimated from seismology and from the theory of the 
figure of the Earth, and (3) that iron would probably be liquid at the temperatures 
and pressures of the core while silicate materials would be solid. 

The hypothesis of an iron core has not gone unchallenged. Kuhn & Rittman 
(1941) propose a core made up of undifferentiated solar matter, arguing that an 
iron core is inconsistent with the relative abundances of the elements. Ramsey 
(1948, 1949, 1950) suggests that the seismic discontinuity at a depth of 2900 km 
represents a phase transition from (Mg, Fe)2SiO4 in the form of olivine in_the 

* Publication No. 108, Institute of Geophysics, University of California. 
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mantle to (Mg, Fe)sSiO, in a metallic phase in the core. Birch (1952) critically 
reviews these suggestions and concludes that the core is principally iron-nickel 
with minor amounts of silicon and carbon since the apparent density of the core 
is somewhat less than the expected density of pure iron. The view that the 2 gookm 
discontinuity is a phase change has many attractive features. However a phase 
change in a multicomponent system almost certainly will be spread out over a 
range of pressures. From the sharpness of the 2g00km discontinuity, a strong 
argument can be made for a marked change in composition at this depth. 

Two classes of observations contain evidence regarding the chemical composi- 
tion of the Earth. Seismic observations and to a lesser degree magnetic and tidal 
observations indicate the physical characteristics of the material making up the 
Earth. The interpretation of these observations in terms of chemical composition 
depends on the measurement or calculation of equations of state of elements and 
compounds valid for temperatures up to a few thousand degrees and pressures 
up to about three million atmospheres. The pressures in the core range from one 
to three million atmospheres. Equations of state for this pressure region have not 
been developed; it is customary to resort to interpolation between observations of 
compression in the laboratory on the one hand and calculation of states of matter 
at pressures so extreme as to cause complex atomic structures to be degenerate on 
the other. 

Direct chemical observations on the outer layer of the Earth, on meteorites 
and on the Sun lead to estimates of the chemical composition of the solar system. 
These estimates can be used to determine the chemical composition of the whole 
Earth, provided some assumption is made regarding the nature of the chemical 
fractionation that has taken place in the solar system. The purpose of the present 
paper is to compare the chemical composition of the core obtained by assuming 
the Earth has the composition of chondritic meteorites with the composition indi- 
cated by seismic observations interpreted in terms of an equation of state inter- 


polated between laboratory data and the solution to the Thomas—Fermi-Dirac 
equation. 


2. Chondritic meteorite model for the Earth 


Urey (1952) and his co-workers (Urey & Craig 1953, Suess & Urey 1956) 
suggest that chondritic meteorites provide a sample of the chemical composition 
of the solar system. Hamaguchi, Reed & Turkevich (1957) have determined the 
uranium content of chondrites using neutron activation techniques. Bate, Huizenga 
& Potratz (1957) have carried out similar studies on thorium. This work combined 
with Edwards’s (1955) determination of the potassium content leads to an estimate 
of the rate at which heat is produced in chondritic matter. 

If the whole Earth is made up of chondritic matter, the heat production should 
be about 2:2 x1029cal/yr. If the Earth contains a chondritic mantle and a non- 
radioactive iron core the heat production should be 1-5 x102%cal/yr (see 
Table 1). 

If the average value of the heat flow from the Earth is taken to be 1-2 x 10-6cm=2s-1 
the rate of heat production in the Earth today is 1-9 x 102° cal/yr, on the assump- 
tion that the Earth is in a thermal steady state. An estimate of the uncertainty in 
this figure comes from two sources: (1) the estimated value of the heat flow may 
be in error and (2) the Earth is not in a thermal steady state. Birch (1954) points o t 
that the heat flow measurements fall within the range (1-2 + 0-4) x 10~®cm~-2s-1 
where we have written the standard error of the measurements listed. The 
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possible uncertainty in the rate of heat production from this source of inaccuracy 
is thus +0-6 x 107° cal/yr. 

If the diffusivity of the mantle is taken to be 0-03 cgs (Clark 1957), then the time 
constant for the mantle of the Earth is much greater than the present age of the 
Earth. We are probably obtaining significant heat flow from only the upper 600 km 
of the mantle. Thus if the radioactivity is uniformly distributed through the mantle 
the Earth’s heat production is actually much greater than that computed from the 


Table 1 
Heat Production 


U X107° g/g 
Th g/g 
K 8-0 x 107* g/g 


Chondritic Mantle and Iron Core 
U 3°29 X 10!* cal/yr 
Th 
K 8-86 x 101° Total 1-53 x 107° cal/yr 


Chondritic Earth 
U 4°82 x 107° cal/yr 
Th 4°56 
K 1°30 X 107° Total 2:23 x 107° cal/yr 


Observed heat flow (1-9+1°0) X 107° cal/yr. 


observed heat flow. Birch (1958) shows that the chondritic model requires the 
concentration of radioactivity in the upper few hundred kilometres of the Earth. 
In order to avoid melting of the mantle below 400 km, MacDonald (1958) suggests 
that the mantle above this depth has the preponderance of the radioactivity of the 
Earth. But this is precisely the region in a more or less thermal steady state with 
the surface. Hence we conclude the uncertainty for the heat production from 
this cause cannot be more than about 0-4 x 1029 cal/yr. 

A problem of somewhat greater difficulty is the uncertainty arising from the 
initial heat distribution in the Earth at the time of its formation. This is of course 
connected with the problem of the present heating or cooling of the mantle. If 
we accept the view that the accreting Earth formed cold and has been heating ever 
since, then we do not have to include this possibility in our calculations of heat 
source production. On the other hand the opposite alternative of origin may give 
us a very large term. It does not seem likely that this figure can be greater than 
0°5 x1079 cal/yr (Jeffreys 1952). 

Combining the several uncertainties, the rate of heat production in the Earth 
can be set at (1°99 + 1-0) x102%cal/yr. The uncertainty in the present rate of heat 
production in the Earth does not allow a distinction between the models of a 
chondritic mantle and a chondritic Earth. The observed heat flow is consistent 
with the assumption that the Earth as a whole is of chondritic composition (Hurley 
1957, Birch 1958). 

Table 2 lists the average composition of chondrites in terms of the number of 
atoms per 108 atoms of silicon. In preparing this table the analysis of Wiik (1956) 
has been added to those reviewed by Urey and Craig. We neglect Urey and Craig’s 


286 

4 

‘ 

‘ 

4 


On the chemical composition of the outer core 287 


separation of chondrites into high and low iron groups. The abundances of ele- 
ments in the Sun’s atmosphere as determined by the Michigan group (Aller 1958) 
are also listed in Table 2 for comparison. 


Table 2 


Chondritic and solar abundances of metallic elements 
in number of atoms per 108 Si atoms 


Element Chondrites Solar atmosphere 


Na 4°94 X 104 5°03 
Mg 9°34 X10 4°80 x 10° 
Al 7°91 X104 4°07 
Si 1-00 X 10° 1:00 X 10° 
4°55 X10° 6-9 
Ss 1:04 108 37 
K 5°94 X 103 2°29 X 10° 
Ca 5°52 X 104 6-03 x 104 
Ti 2°27 X 108 2°29 X 108 
Cr 7°70 X10 2°51 
Mn 5°64 x 10° 5°30 X 10° 
Fe 7°12 X10 1°45 X10° 
Co 2°21 X70" 
Ni 3°64 x 104 1°59 X 104 


The data of Table 2 can be used to illustrate the inconsistency of assuming 
a mantle made up largely of olivine and a core of iron and assuming a chondritic 
Earth. The three principal metals in chondrites and in the Sun are silicon, iron, 
and magnesium. If we assume a mantle made up of olivine (Mgo.9Feo.1)2SiO4 and 
a core of iron-nickel, the ratio of the sum of the iron, magnesium and nickel 
abundances to that of silicon is 3-2 as compared to the chondritic ratio of 1°65 
(Table 3). A similar contradiction arises if we assume the mantle is a peridotite. 


Table 3 
Atomic ratios of iron, magnesium and silicon in various earth models 
Model (Fe +Mg +Ni)/Si Fe/Mg 
Composition Composition for whole for whole 
of Mantle of Core Earth Earth 
(Mg o- 9Fe o-1)25iOz (Fe,5Ni) 3°2 0-78 
9Feo.1)SiOs (Fe,;Ni) 2°2 1-04 
Peridotite (Fe,;Ni) 2°9 
Chondritic meteorites 1°65 o-71 
Solar atmosphere o-71 0°30 


If we assume the mantle is made up of peridotite having the composition given 

by Daly (1942) (Table 3) and we assume an iron-nickel core, the abundance ratio 

(Fe+ Ni+Mg)/Si for the earth as a whole is 2-9. Maintaining the hypothesis 

that chondrites are representative chemical samples of the solar system we are 

left with several possibilities: (1) the composition of the mantle, or core, or both 
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assumed above is incorrect; (2) some process of differentiation may have concen- 
trated iron and magnesium relative to silicon during the Earth’s formation, at the 
same time keeping the concentration of radioactive elements constant. Ruling out 
the second hypothesis as too restrictive, we examine alternative compositions for 
the core and mantle. 

An iron core is consistent with a chondritic composition for the Earth, provided 
the ratio (Fe+ Mg)/Si in the mantle is about unity, as illustrated in Table 3. This 
implies a mantle dominantly of pyroxene composition. The seismic velocities 
observed in the upper mantle rule out enstatite (MgSiOg) or hypersthene 
(Mgo.9Feo.1)SiOg as major constituents since the quantity (0P/0p), = ¢ for these 
minerals (Table 4) is much lower than that observed at depths of 200 to 400 km 


Table 4 
Values of 

(Kr/P)o (Ks/P) 
Enstatite MgSiO, 29 
Hypersthene (Mg,Fe)SiO; 28 
Forsterite Mg,SiO, 36 
100 km depth in Earth 36°8 
200 km 40°1 
300 km 43°2 
400 km 47°0 


(Birch 1952). It is possible that pyroxenes undergo a phase transition to a form 
with higher ¢ at pressures less than those corresponding to depths of 100-200 km 
but this transition has not been observed experimentally (MacDonald 1957). A 
mixture of iron-poor pyroxenes and garnets (eclogite) may give a suitable value of 
¢ for the upper mantle (Birch 1952). 

On the basis of the solid nature of the upper mantle, MacDonald (1958) 
shows that the radioactive elements must be concentrated in, at most, the upper 
400km of the mantle. Assuming that Na and Al are concentrated along with 
K, U, and Th in this region the material of the upper 400 km of the mantle would 
have the composition of an olivine gabbro or eclogite with seismic velocities requir- 
ing a mineral assemblage of olivine, magnesium-rich garnet and alkali-rich 
pyroxene. Assuming an upper mantle (400km) of eclogite of composition given 
in Table 5 and a lower mantle of dunite, an iron core leads to a (Fe + Mg)/Si ratio 
for the Earth of 2-9, and Fe/Mg ratio of 0-94, both values in disagreement with the 
chondritic ratios. 

A chondritic model containing an upper mantle of eclogite and a lower mantle 
of dunite requires a core having a composition (Fe, Ni);.gSi. The iron-magnesium 
ratio in such a model is Fe/Mg = 0-71. The ratio of iron to silicon in the core is 
raised, of course, if a higher silicon content is assumed for the mantle. The com- 
position of the core is not uniquely determined by abundance arguments because 
of the uncertain composition of the mantle. We do conclude that a core composi- 
tion of (Fe, Ni);.¢_2.0Si is consistent with the chondritic model for the Earth 
and with a eclogite-dunite model for the mantle. If a major part of the mantle 
has the composition of eclogite with (Fe+Mg)/Si~ 1, the amount of silicon 
required for the core is reduced in proportion to the ratio of pyroxene-garnet 
to olivine. Such a model for the mantle is inconsistent with a pure iron core as 
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is any model for the mantle in which the (Fe + Mg)/Si ratio in the mantle is appre- 
ciably greater than one. An iron-nickel-silicon alloy cannot be ruled out on chemi- 
cal relations observed at one atmosphere, since iron and silicon are completely 
miscible in the liquid state and even in the solid phase the solubility of silicon in 
iron ranges from Feg3.;Si at 1 000 °C to Feg.gSi at 400 °C (Greines, Marsh & Staugh- 
ton 1933; Farquhar, Lipson & Weill 1945). 


Table 5 
Composition of possible components of mantle in number of atoms per 108 atoms of Si 

eclogite! peridotite?* dunite?* 
Na 5°82 x 104 2°78 x 104 2°39 X 10° 
Mg 4°44 x10° 1°25 X10° 1°70 X10 
Al 2°54 1°29 X 10° 2°50 
Si 1-00 X 10° 1-00 X 10° 1°00 X 10% 
1°46 x 10° 1°92 X 10° X10? 
K 1°49 X 104 7°10 X 103 6:29 10? 
Ca 2:21 X 108 8-70 x 104 1°85 x 104 
Ti 1°36 x 104 1°71 X 103 3°71 X10? 
Cr 1°36 x 108 
Mn 3'29 X10 3°66 x 10% 3°34 X10 
Fe 2°41 X 10° 1°58 x10° 1°67 x 10° 
Ni 6°94 X 10? 


1 Average of 11 non-kyanite bearing eclogites found as inclusions in kimberlite pipes in 
South Africa (Williams 1932). 
2 After Daly (1942). 


3- Calculation of the representative atomic number 


Various attempts have been made to apply the Thomas—Fermi equation of state 
to elements at conditions prevailing in the core (Jensen 1938, Elsasser 1951, 
Knopoff & Uffen 1954). Marshak & Bethe (1940) indicate the solution to the Tho- 
mas—Fermi equation for compounds. Knopoff and Uffen show that the solution 
to the Thomas—Fermi-Dirac (T.F.D.) equation for compounds simplifies if the 
compound is characterized by a single “representative atomic number” (Z) which 
is a weighted mean atomic number for the compound. The equation of state for 
the compound is the T.F.D. equation of a pure element having the Z of the com- 
pound. No difficulty results if the Z assumes non-integer values. The Z is 
calculated at pressures at which the T.F.D. equation applies but the number is 
associated with the compound in the intermediate pressure range. We will use 
the T.F.D. equation to obtain the expression for Z in terms of the atomic number 
of the constituent elements. The actual equation of state applied to the core is an 
interpolation between the Birch-Murnaghan equation applicable at pressures up 
to at least 0-5 x10%atm (Duvall & Zwolinsky 1955) and the T.F.D. equation 
applicable at 107 atm and higher. 

Knopoff and Uffen calculated Z by numerically comparing equations of state 
of assemblages of elements with the equation of state of pure elements. March 
(1955) and Gilvarry (1957) present asymptotic solutions to the T.F.D. equation 
which suggest that the calculation of Z can proceed analytically rather than semi- 
empirically, provided the calculation is carried out at pressures sufficiently high 
for the asymptotic solution to hold. These pressures are somewhat greater than 
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the lowest pressures at which the T.F.D. model is assumed valid, that is to say 
107atm. The asymptotic formula accurately represents the equation of state of 
elements at pressures of 109 atm and greater. 

Neglecting binding energies and assuming the equality of mass of proton and 
neutron, the equation of state valid at extreme pressures and absolute zero is 


P = pressure 

Z = atomic number 

M = atomic weight 

= density 

= 
= 

= 

= Planck’s constant 

= electronic charge 

= mass of the electron 
= ratio of mass of proton to mass of electron. 


where 


We consider a compound Y with a chemical formula 
Xn} Xn?...Xn™ 


where the n; are the number of atoms of the element X‘ in the chemical formula. 
The condition that a compound with representative atomic number Z and repre- 
sentative atomic weight M has the same density as Y is then 


+ + = 


I 

where Z; is the atomic number of element X“), M;, is the atomic weight of element 

We now make the assumption that the abundance of neutrons in our hypo- 
thetical compound with representative atomic number Z is determined by the 
actual abundance of neutrons in compound Y. We state that the ratio of neutrons 
to protons in both Y and in the equivalent compound is the same. 


™ 

mM 
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Combining equations (1) and (2) we obtain an expression for Z, viz: 


2; 

ia 
223 = 


m 
Z, 


" 

4 


On the chemical composition of the outer core 291 


Representative values of Z for various materials are listed in Table 6. The 
impressive feature of the values listed in Table 6 is.that Z for all silicates and for 
all oxides except those of iron is well below 20. 


Table 6 
Representative atomic numbers of compounds 
Material Formula 

Quartz SiO, 10.7 
Anorthite CaAl,Si,0. 
Orthoclase KAISi;0, 11.6 
Albite NaAlSi;0, 10.6 
Diopside CaMgSi.O, 12.0 
Enstatite MgSiO; 10°5 
Jadeite NaAlSi,O, 10°5 
Fayalite Fe,SiO, 17°7 
Forsterite Mg,SiO, 10°5 
Andradite CasFe2(SiO,)s 
Grossularite CasAl,(SiO,)s 12°6 
Almandite FesAl3(SiO,)s 14°8 
Pyrope Mg:Al,(SiO,): 11°7 
Spodumene LiAlSi,O, 
Beryl 
Corundum Al,O; 10°4 
Wustite FeO 21'2 
Magnetite Fe;0, 20°1 
Hematite Fe,0; 19°6 
Periclase MgO 10°3 
Spinel MgAl,O, 10'5 
Rutile TiO, 15°5 
Oldhamite CaS 18-2 
Troilite FeS 22:0 
Pyrite FeS; 20°3 
Galena PbS 69:0 
Sphalerite ZnS 24°8 
Iron Fe 26 
Nickel Ni 28 

Fe,Si 

Fei..Si 


4- Chemical composition of the core derived from seismology and equa- 
tion of state arguments 


On the basis of Bullen’s density distribution, various workers infer that the 
outer core is somewhat less dense than pure iron or an iron-nickel alloy (Birch 
1952, Bullen 1952). Knopoff and Uffen, using an interpolated equation of state, 
obtain a Z for the outer core of 22. The exact value of Z is uncertain because of 
the nature of the interpolation and the small but unknown effect of temperature 
on the equation of state. It appears unlikely that the uncertainty could be as great 
as four units in Z but might be as large as one or two units. We will use a value of 
22 for Z for the outer core. Since the atomic numbers of iron and nickel are 26 
and 28 respectively, some lighter elements must be present in the core in order to 
bring Z down to 22. 

The abundance arguments reviewed above suggest that silicon might be a 
major constituent of the Earth’s core. In Figure 1 we show the dependence of Z 
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for a solution of iron and silicon on the amount of silicon present. If Z for the core 
lies between 22 and 23, then the outer core must contain between 20 and 30 
weight per cent silicon corresponding to the formulae Fej.2Si(Z = 22) and Fep.;Si 
(Z = 23) respectively. The composition of Fe;.gSi suggested by the abundance 
of elements is thus in excellent agreement with the composition deduced on the 
basis of equations of state, provided silicon is assumed to be the alloying element. 
The equation of state cannot rule out the possibility of a liquid solution of 
Fe-Si-Mg. Figure 2 gives the relations between the number of atoms of these 
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Fic. 1.—The number per cent of atoms and the weight per cent of elements 


required to dilute iron to a given value of Z for the systems Fe,Hy, Fe,Hey, 
Fe,Siy, Fe,S,. 
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three constituents and Z. Since silicon and magnesium have atomic numbers which 
are close to one another, we see that the range Fej.9(Si, Mg) to Fee.;(Si, Mg) fits 
estimates of Z. Similarly nickel and iron have similar atomic numbers and solu- 
tions of the form (Fe, Ni):.2Si to (Fe, Ni)z.1Si fit the requirements of the equation 
of state. The abundance of Ni in chondrites suggests a Fe/Ni ratio of about 20. 

The limits imposed by the estimated values of Z do not uniquely determine the 
composition of the core. Other elements and combination of elements can be used 
to reduce the atomic number of the material in the core as is illustrated in Figure 1. 
Solar and chondritic abundances suggest that hydrogen, helium, carbon, oxygen, 
nitrogen, magnesium and sulphur are likely candidates for reducing the value 
of Z of the core. All but magnesium are either highly volatile or form highly volatile 
compounds under a wide range of conditions. If we consider chondritic meteorites 
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as indicating those elements other than silicon that might be retained by a planet 
in large quantities then only oxygen, magnesium and sulphur are possible diluents 
of iron-nickel. 

Any high-pressure modification of a silicate material is ruled out as the sole con- 
stituent of the outer core, since the Z values for silicates are all considerably lower 
than 22. A mixture of magnesium silicate and metallic iron is a possibility. If 
such a silicate were initially dispersed in a molten iron core it would almost 
certainly gravitate towards the outer edge of the core and merge with material 
of the mantle. Alternatively, the silicate material might be liquid. In this case we 
would have to assume that the metal phase and silicate phase were miscible, other- 
wise these two phases would tend to separate under gravity into a region of high Z 


Si (Z= 14) o 


Fic. 2.—The fraction of number of atoms for different values of Z in the 
system Fe,Mg,Si,. 


below a region at the top of the core of lower Z. Such a liquid separation in the 
outer core contradicts the seismic evidence. Molten silicates and metallic iron are 
immiscible at surface pressures; there is no reason to expect a marked increase of 
solubility at core conditions. 

As is evident from Figure 1, a large number of atoms but only a small weight 
of either hydrogen or helium reduces Z for the core to the range of 22-23. Both 
of these elements are abundant in the Sun but are not major constituents of 
meteorites. There is no guarantee that the present gaseous composition of 
meteorites in any way represents their original volatile element composition so 
that hydrogen and helium remain possible if perhaps unlikely diluents of the 
core. 

Sulphur is a possible major component of the core since FeS has a Z of 22. 
Iron sulphide is found in meteorites as troilite. The abundance of sulphur in 
chondrites is about a tenth that of silicon. If the core were FeS and the mantle 
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were eclogite—dunite with no sulphides, then the ratio of sulphur to silicon would be 
S/Si = 0-77 rather than o-1 or 0-37 as indicated by chondritic and solar abundances 
respectively. A postulate of an iron-sulphur mixture with the appropriate Z for 
the core requires an initial high abundance of sulphur or some process by which 
silicon is lost to the Earth with respect to sulphur. Neither of these hypotheses 
seems probable though they cannot i. ruled out completely. On the chondritic 
model of the Earth it is possible to account for the sulphur by 3 to 4 per cent by 
weight of metallic sulphides in the mantle with the remaining sulphur in the core. 
However sulphur cannot be ruled out as a possible component of the core though 


quantitatively the abundance of available sulphur is too low for it to be the major 
diluent of the core. 


5. Problems raised by a core of iron silicide 


We have shown that a core of composition (Fe, Ni);.gSi is consistent with a 
chondritic composition for the Earth, provided the mantle is made up of an upper 
layer having the composition of eclogite and a lower shell of dunite composition. 
The Z of Fej.gSi is 22°8 while that calculated for the core from Bullen’s density 
data is 22. The agreement between two methods of estimating the composition of 
the core may be fortuitous. Since the agreement is sufficiently striking we examine 
two of the problems raised by the postulate of a core containing abundant silicon. 

The dynamo theories of the Earth’s magnetic field require a material of rela- 
tively high electrical conductivity in the core. Runcorn (1956) on the basis of the 
form of the secular variation spectrum derives a value of 1oohm~!cm~! as the 
conductivity of the core. The analysis by Bullard, Freedman, Gellman & Nixon 
(1950) of the westward drift shows that a value of 3 x 10 ohm~! cm~? is a reasonable 
estimate of the conductivity and it is this latter value that is generally used in dis- 
cussions of the theory of the magnetic field. Various attempts have been made to 
estimate the conductivity of iron at conditions prevailing in the core. Elsasser’s 
(1950) estimate is 1-3 x104ohm=!cm7! for iron at gooo °C or 3 x104 at a more 
reasonable temperature of 4000°C. A conductivity as high as 3 x 104ohm=!cm7 
requires that the currents producing the time-dependent part of the field originate 
in a very thin layer near the mantle boundary. The electrical and magnetic 
properties of iron-silicon alloys at surface conditions are well known because of 
the use of iron-silicon alloys as transformer cores and in other electrical equipment. 
The conductivity of pure iron (1 x10°ohm=!cm-~!) decreases markedly with the 
addition of silicon to the value of 1 x 104 ohm=!cm-~ at 11 per cent silicon by weight. 
At 14°5 per cent silicon the conductivity reaches a maximum of 3 x 104 ohm~! cm7! 
and decreases with further addition of silicon (Bozorth 1951). The extrapolation 
of these conductivity values to conditions of the core is most uncertain but it 
seems likely that the conductivity of Fe;.gSi would be lower than that of pure iron 
by a factor of three and perhaps lower by an order of magnitude. Both theory 
and “observations” require a conductivity somewhat less than that estimated for 
pure iron so that a core composition of Fe;.gSi is consistent with what is presently 
known about the Earth’s magnetic field. 

The existence of iron meteorites free of iron silicide suggests a further problem. 
The hypothesis of an iron-silicon core requires that the process by which iron is 
freed from oxygen in the Earth differ from the mechanism by which free iron is 
produced in meteorites. The closeness of the ratio of iron to silicon in the postulated 
core and in the mineral fayalite provides the tantalizing suggestion that the core 
is produced by the squeezing of oxygen out of a silicate having a composition 
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Fe2SiO4. A requirement for such a process is that the system be open to oxygen 
so that the volume change in the reaction FegSiO4 = FegSi+ 2Os is solely the differ- 
ence in volume of FegSi and Fe2SiO4. A convective process sweeping away the 
volatile components might provide the mechanism by which the system is opened 
to oxygen. At zero oxygen pressure FegSi+2Qz is stable relative to FegSiO4 as 
measured by the heat of formation of Fe2SiO4. Neglecting the changes of entropy 
and volume with temperature and pressure, the change in the Gibbs free energy in 
the reaction 


FeeSiO4 = FeeSi+ 202 
at high total pressure is 
AG = AGo+ PAV;—(T—T)AS 


where AGpo is the change in the Gibbs free energy for the reaction at surface 
conditions, AV, is the difference in the volumes of the solid phases and AS is the 
total change of entropy. Assuming a partial pressure of oxygen of 1 atm and standard 
thermochemical data (Kubaschewski & Evans 1951) the pressure needed to con- 
vert FegSiO4 into Fe2Si plus oxygen is of the order of 7 x 105atm while 5 x 106 atm 
are needed to squeeze oxygen out of forsterite under the same conditions. These 
calculations neglect the effect of the possible phase changes of olivine to spinel 
but since the surface oxygen pressure is of the order of 1 atm they indicate the 
minimum pressures required to break down iron silicates. These pressures are 
reached within the Earth but probably not in the planet from which the meteorites 
may have been derived. 

The suggestion that the outer core be a silicide has also been made by Ring- 
wood (in press). Ringwood has approached the problem from the point of view that 
the solar matter undergoes progressively greater chemical reduction in depth in 
the Earth. Although the general conclusions are similar the process of introducing 
a silicide into the core is different; Ringwood suggests an amount of silicon similar 
to that proposed here. 


Institute of Geophysics, 
University of California, 
Los Angeles, 
Cal., U.S.A. 
1958 May 23 
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A New Determination of the Differences in Gravity 
between the N.P.L. Teddington and the South African 
Fundamental and Secondary Stations 


D. I. Gough* 


(Received 1958 July 10) 


Summary 


Observations have been made with the Cambridge pendulum 
apparatus at the British Fundamental gravity station at the N.P.L. 
Teddington, at the South African Fundamental gravity station at 
Johannesburg, and at the South African secondary stations at Pretoria 
and Mowbray. Two pendulums were exposed to mercury during the 
project. The remaining four, namely IA, IB, IC and VIA, give 
for the mean gravity differences 


Teddington—-Mowbray 1°5495 + 0-0001 cm/s? 
Teddington—Johannesburg = 2°6472 + 0:0003 cm/s? 
Teddington-Pretoria = 2°5672 + 0:0003 cm/s? 
whence, adopting g = 981-1963 cm/s? at Teddington, 
g& = 979°6468 cm/s? at Mowbray 
978-5491 at Johannesburg 
978-6291 at Pretoria. 


These values agree well with the determinations of 1948-9 by Hales and 
Gough. 


1. Introduction 


The gravity survey of the Union of South Africa has recently been completed 
with the publication by the Union Geological Survey of a Bouguer anomaly mapt 
overprinted on the 1 : 1 M geological map of the Union. The gravity map repre- 
sents gravity meter observations at about 6 300 stations. The gravity meter observa- 
tions were controlled by the primary network of pendulum stations, forty in the 
Union and thirteen in neighbouring territories, observed by Hales and Gough in 
1948-9 using the Cambridge invar pendulums (Hales & Gough 1950). The pendu- 
lum stations in Southern Africa were occupied in loops beginning and ending at 
the Fundamental station at the Bernard Price Institute of Geophysical Research, 
Johannesburg, or at the Trigonometrical Survey Office, Mowbray, Cape Town. 
In consequence, measurements were made several times during the pendulum 


* National Physical Research Laboratory, Council for Scientific and Industrial Research, 
Pretoria. 

t No formal decision has been taken that Johannesburg shall be the South African Fundamental 
Station, but it is the South African first order station in the international network. 

~ This map may be purchased from the Government Printer, Pretoria. 
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survey at each of these stations, and the difference in gravity between Cambridge 
and each of them was well determined. A re-connexién between Mowbray and 
Johannesburg and a first-order station in another continent is therefore equivalent 
to a re-connexion between the national survey and the international gravity 
network. 

In addition to the general desirability of a new determination of the differences 
between the British Fundamental station at the National Physical Laboratory, 
Teddington, and the South African stations, there was a specific reason for a new 
connexion. This was the uncertainty regarding the effects of magnetic interaction 
between the invar pendulums and the geomagnetic field, or between the pendulums 
and the mu-metal shield with which the pendulum swinging case was lined 
during the 1948-9 work (Bullard 1933, Cook 1954). This uncertainty was increased 
by the reversal of the vertical component of the geomagnetic field between Britain 
and South Africa. Further, the Cambridge pendulums have on various occasions 
shown small discrepancies in the measured gravity differences given by Set I 
and Set VI respectively (Cook 1954) and this discrepancy was exceptionally large 
for the 1948-9 survey (Hales & Gough 1950, Cook 1957). Cook has advanced 
an explanation of this discrepancy based on magnetic perturbation (Cook 
1954): 

For these reasons the writer accepted with pleasure an offer by Mr B. C. 
Browne to make the Cambridge pendulum apparatus available for a new connexion 
between Teddington and the South African stations. Mr J. E. Jackson and the 
writer made the initial observations at Teddington and the apparatus travelled 
to Cape Town in the writer’s care. After the South African observations had been 
completed, the apparatus was shipped to Southampton and Mr Browne made the 
closing observations at Teddington. 


2, Measurement procedure 


The positions of the South African stations are given in Table 1. It will be 
seen that measurements were made at the Transvaal Museum, Pretoria, station 


Table 1 


Positions of stations in S. Africa 


Designation Description Latitude Longitude Height 
s E (m) 

Mowbray Trigonometrical Survey 33°57" 18°28’+1 38°4 
Office, Mowbray, 
Cape Town 

Johannesburg Bernard Price Institute, 26°11''5 28°o1’.8 1755°0 
Johannesburg. 

Pretoria Transvaal Museum, 25°45 28°11'*4 1337°1 
Pretoria 


as well as at Mowbray and Johannesburg. These were made because the gravity 
meter network of the gravity survey used Pretoria as origin. The stations were 
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observed in the sequence Teddington-Mowbray-Johannesburg-—Pretoria (two sets 
of 

The Cambridge pendulum apparatus has been described frequently (see 
Browne, Cook, McCarthy & Parasnis 1950 for references). Two invar pendulums 
of half-period slightly greater than 0-5s were swung in antiphase about an axis 
which was along the magnetic meridian. The mu-metal lining formerly used to 
reduce the field within the swinging case had been removed, and instead a Helm- 
holtz pair of coils annulled the vertical component of the geomagnetic field, 
The pendulum periods were measured by comparison with a clock driven by a 
quartz crystal oscillator of frequency about 100000¢/s, light signals from the 
pendulums and clock being recorded together on a chronograph drum. The 
air pressure in the swinging case varied from 18-9 to 27-2 mm of mercury, and was 
measured throughout by one U-tube manometer inside the swinging case, checked 
through an external connexion against the N.P.L. standard U-tube at Teddington, 
and against a second manometer in South Africa. At Teddington the quartz oscil- 
lator was rated by beating the second harmonic of its output with the 200 kc/s 
carrier of the B.B.C. Light Programme. For rating the clock in South Africa, Mr 
L. Flavell of the Department of Geodesy and Geophysics, Cambridge, attached to 
the clock a cam-operated contact-breaker which enabled the output of a radio 
receiver to be interrupted once every half-second of clock time. The contact- 
breaker could be rotated round the axis of the cam and its angular displacement 
was shown on a dial in hundredths of a half-second of clock time. The device 
was used with the 5 Mc/s transmission of ZUO, the radio station of the Union 
Observatory, Johannesburg. This transmission included “ticks” every second 
consisting of 10 cycles of 1o000c/s tone. The contact-breaker was turned until 
the “ticks” just disappeared, the dial was read, and the time was noted. By 
repeating the setting several times over several hours, the rate of the clock 
could be found with a precision limited by radio interference, but always in 
the present work high enough to make error in clock rate negligible compared 
with the overall variations in pendulum periods. In practice, the contact- 
breaker was set so that the “tick” just disappeared inside the leading and 
trailing edges of the silent interval, and the mean of the two dial settings was 
used, so as to minimize effects of variation in signal to noise ratio and of 
eccentricity of the contact-breaker. Several such pairs of settings made up an 
observation. 

Both at Teddington and in South Africa rating observations were made at 
intervals during each day on which pendulums were swung. At Mowbray, radio 
reception conditions made it necessary to work at night. The rates measured 
are given in Table 2. Most of the variation in rate resulted from temperature 
variation. 

One set of observations at a station consisted of three swings (in two cases 
four swings) of each of four pairs of pendulums, the pairs being I A, B; I A, C; 
VIA, B and VIB, C. In most cases only one pair of pendulums was swung in a 
day, to give maximum time for settling down of temperature after handling. 
One swing was about 7000 half-seconds long, all other swings had duration 
about 8 ooo half-seconds. 

Before the first observations at Mowbray an externally connected manometer 
broke in such a way that mercury was forced into the swinging case by the ingoing 
air and sprayed over the VIB and VIC pendulums. The mercury droplets did 
not appear to wet the gold plating, and when the pendulums had been wiped they 
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Table 2 
Rate of quartz oscillator © 

Pend- Rate Pend- Rate 
Date ulums Date ulums x 108 

in use in use 

Teddington Pretoria 
20/9/56 VI BC +529+1 17/11/56 I AB +480+2 
21 VI AB +531+2 19 I AC +480+4 
21 I AB +530+1 20 VI AB +47043 
22 I AB +532+1 21 VI BC +469+3 
24 I AC +s525+1 
Mowbray Johannesburg — 
19-20/10/56 I AB +517+43 23/11/56 I AB +514+43 
20-21 I AC +509+6 26 I AC +510+2 
23-24 VI BC +523+5 27 VI BC +508+2 
26-27 VI AB +535+5 28 VI AB +506+2 
Johannesburg Mowbray 
§-6/11/56 VI BC +514+1 1-2/12/56 VI BC +498+3 
7 VI AB +525+3 2-3 VI AB +498+3 
8 I AB +520+2 3-4 I AB +495+5 
9 I AC +502+4 455 I AC +s501+7 
Pretoria Teddington 

13/11/56 I AC +503+3 20/2/57 I AB +520+4 
14 I AB +489+4 21 I AC +518 
15 VI BC +494+2 21 VI AB +515+5 
16 VI AB +456+3 22 VI BC +517+6 


Note—The errors quoted are mean square errors. 


appeared unchanged. During subsequent weeks, however, spots and cloudy 
areas of amalgam appeared, especially on pendulum VIB. 

The usual corrections were applied to yield half-periods at zero amplitude and 
pressure, o°C and zero clock rate. The amplitude correction varied from — 143 
to — 174, the pressure correction from —279 to — 369, the temperature correction 
from —468 to —629 and the clock rate corrections from —231 to —271, all in 
units of 10-8s as corrections to the pendulum half-period. 


3. The observations 


The corrected pendulum half-periods are given in Table 3. It will be seen that 
there appears to be a general tendency for the half-periods of all the pendulums 
except VI B to increase between the first and second measurements at each station. 
In particular at Pretoria the mean period for the second set of observations is 
greater than that for the first set for all six pendulums, the differences in half- 
period ranging from +9 to +32 x 10-8s. The mean difference is 20 x 1078s, 
corresponding to 0-8 mgal. Between these two sets of observations, which together 
occupied nine days, the pendulum swinging case was lifted, re-oriented with 
respect to the magnetic meridian and re-levelled. The pendulums were handled 
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with the usual stringent precautions, and in any case it is unlikely that all six in- 
creased in length. The physical conditions associated with corrections changed only 
slightly. The magnetic moments of the pendulums changed only slightly. Possible 
levelling errors were far too small to account for the differences. The data suggest 
that there are unknown effects associated with setting up the swinging case which 
may cause variations of the order of one milligal. If so, the variation of period 
between swings of a set with the swinging case undisturbed will not reflect the 
true precision of a measurement. It might be useful, though tedious, to study 
the effect by means of a series of sets of swings at some one station. 

The closure differences in half-period at Teddington range from +13 to +56 
x 1078s for the four pendulums used in deriving the differences between Tedding- 
ton and the South African stations. The mean closure difference of +30 x 1078s 
corresponds to 1-2 mgal. 

The VI pendulums were swung in pairs VIA, B and VIB, C instead of 
VIA, B and VIA, C, on the suggestion of Mr Browne, so that the VI B pendulum 
was used on both left and right pairs of agates, and any change in the difference 
between the rigidities of the two supports would show as a change in the difference 
in the half-periods. Part of the difference will be due to sway coupling to VIA 
and VIC respectively, but this will be constant (Browne & others 1950). Table 4 
shows that there was such a difference and that it was apparently not constant. 


Table 4 


Difference, A, of mean half-period of Pendulum VIB on right-hand agate from that 
on left-hand agate 


Station M J 4 J M 


A(10-8s) 36 52 31 27 23 21 19 II 


The gravity differences from Teddington are given in Table 5. These have 
been calculated using the equations 


Ag = —3°87634 x10 5AT+ 1-149 x 10712(AT)? 
for the set I and 


Ag = —3°88237 x10 5AT+ 1-152 


for the set VI. Gravity differences are given for each of the six pendulums. Mean 
differences are quoted for all six pendulums, and for LA, 1B, IC and VIA. It 
is the writer’s view that the latter should be adopted in view of the accident with 
the mercury. Any effect of the mercury on the half-periods of pendulums VI B 
and VIC is clearly of the same order as changes which occurred in those of the 
other pendulums, as Table 5 shows. Nevertheless pendulums VIB and VIC 
are known to have suffered a physical change and it is therefore considered better 
to exclude them. All six pendulums were used in deriving the mean differences 
between South African stations, as the mercury accident had occurred before the 
first South African observations. 
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Table 5 
Gravity differences 


Unit cm/s? 

Pendulum T-M T-J T-P M-J P-J 
IA 1°5494 26476 2°5673 1:0982 00803 
IB 1°5494 26468 2°5669 1:0974 0:0799 
IC 1°5495 2°6471 2:5669 1:0976 00802 
VIA 1°5498 26475 2°5678 10977 0:0797 
VIB 1°5485 2°6462 2°5663 1:0977 0:0799 
VIC 1°5487 26466 2°5670 1°0979 0:0796 

Mean of six 1°5492 2°6470 2°5670 1:0978 00799 

Std. error of mean 00003 0:0003 0-0002 00001 00001 

Mean of [A, IB, IC, VIA 1°5495 2°6472 2°5672 

Std. error of mean 0-000! 0:0003 0:0003 


T = Teddington, M = Mowbray, J = Johannesburg, P = Pretoria. 


Browne & others (1950) point out that effects of sway of the support and of 
ground motion are less if gravity differences are calculated from the mean of the 
periods of the two pendulums swung as a pair, than if individual pendulum periods 
are used. This involves giving double weight to the pendulum in each set which 
is swung as a member of both pairs (I A and VIB in the present work). Professor 
Hales and the writer are of the opinion, as at the time of the earlier pendulum 
measurements, that where all observation sites are stable in character, sway and 
ground motion are likely to be less important than changes in length of the pen- 
dulums, and unknown effects such as that discussed in connexion with the observa- 
tions at Pretoria. If this is so then, as Cook has remarked (Cook 1957, p. 232), it 
it best to use individual pendulum periods with unit weight in calculating gravity 
differences. Gravity differences have therefore been found from the individual 
pendulum half-periods with unit weight. The alternative decision would, in view 
of the mercury accident, have made the connexion with Teddington rest upon 
Set I of pendulums alone. 

Table 6 gives values of gravity at the South African stations, taking the value to 
be 981-1963 cm/s? at Teddington (Cook 1952). Values given by Rose & Woollard 


Table 6 
Recent pendulum determinations at the South African Gravity Base Stations 


Unit cm/s? 
Date Reference Base Pendulums Johannes- Mowbray Pretoria 
burg 
1948-9 Hales & Gough Cambridge Cambridge 978-5497 979°6474 978°6293 
1950 invar 
1955 Rose & Woollard Washington Gulf “M” 978-5504  979:6486 
1956 quartz 


1956 Gough (present work) Teddington Cambridge 978-5491 979°6468 978-6291 
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(1956) and by Hales & Gough (1950) are quoted for comparison. The last mentioned 
have been increased by 3-5 mgal to correct for the change in the assumed value of 
gravity to 981-2685 cm/s? in Pendulum House, Cambridge (Cook 1953). It will be 
seen that agreement is satisfactory between the two sets of data from the Cambridge 
pendulums, and that magnetic perturbations in the earlier work do not appear to 
have affected the mean differences greatly. Agreement with the Gulf quartz pen- 
dulums is less close but this may reflect systematic differences between the two sets 
of apparatus. It should be noted that Rose and Woollard state, in the preface to 
the Technical Report quoted, that some of the results cited in this report, particu- 
larly pendulum results in South Africa, are of a preliminary nature and are not to 
be regarded as final values. Rose & Woollard (1956) give results from two con- 
nexions between Washington and South African stations using Worden gravi- 
meters. The values given agree reasonably with those quoted in Table 6. 

Table 7 gives gravity differences within South Africa found in 1948-9 and in 
the present work. The difference between Pretoria and Johannesburg from the 


Table 7 
Gravity differences within South Africa in 1948-9 and in 1956. 
Unit cm/s? 
M-J P-J 
Hales & Gough 1950 1°0977 +0°0003 0:0800 + 0:0008 
Present work 1:0978 +0-0001 0:0799 +0:0001 
= Johannesburg M = Mowbray P = Pretoria 


earlier work is given as 80-0 mgal in Table 7 and as 79-6 mgal by Hales & Gough 
(1950). The latter value was derived from thirty-six swings made over a period 
of four months at Johannesburg and eight swings made on 1948 August 5 and 6 
at Pretoria. A difference more comparable with the difference found in the 
present work is that between the 10 swings at Johannesburg between 1948 August 
3 and g and the 8 swings at Pretoria. This is the difference quoted in Table 7. 
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Magnetic Cleaning of Rocks in Palaeomagnetic Research 
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Summary 


It has been observed that the remanent magnetism of rocks is often 
the resultant of two magnetizations: firstly a magnetization with low 
coercive force, mostly directed by the local present Earth’s field; 
and secondly a magnetization with high coercive force which shows 
the direction of the field from the period of the rock formation. 

It is shown that it is possible to remove the disturbing magnetization 
of low coercive force while saving a good deal of the fossil magnetization 
with the aid of partial demagnetization. This was carried out by a 
combination of alternating fields and heating. 

Generally a demagnetization with alternating fields up to 300 oersted 


(effective) at a temperature of 150°C appeared to be sufficient to clean 
the rocks. 


1. Introduction 


In palaeomagnetic investigations it has often been observed that there is a great 
dispersion in the directions of magnetization of rock-specimens of the same age 
and the same region even if the specimens come from the same formation. In 
this case the directions plotted in a stereographic projection often show a belt and 
are situated between the directions of the ancient field and the Earth’s present 
field in that region. This has been already reported by some investigators and 
the explanation introduced is that the present field has influenced the magnetism 
of the rocks (Green 1957 b). 

To check whether the magnetization in rocks is a stable fossil magnetization, 
some stability tests have already been developed (Graham 1955, Irving & Runcorn 
1957 and Nagata 1953). However, these stability tests give only a general idea on 
the reliability of the magnetism of rocks. In practice it is impossible to perform 
these tests on all rock-specimens of the collection. For this reason an attempt has 
been made to develop a method of magnetic cleaning which, after application, 
leaves only the stable fossil remanent magnetism. 


2. Principle of magnetic cleaning 


Our starting-point was the idea that the magnetization measured in rocks is a 
combination of different magnetizations, namely a magnetization directed according 
to the ancient field in the period of deposition, metamorphism or diagenesis of 
the rocks and magnetizations of later periods. It is known that the magnetic 


* Koninklijk Nederlands Meteorologisch Instituut, De Bilt. 
t Mineralogisch-Geologisch Instituut, Rijksuniversiteit, Utrecht. 
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particles in a stone have different coercive forces dependent on their composition, 
shape and size. The magnetization of the particles with a fairly high coercive 
force can only be changed by a fairly strong magnetic field. In palaeomagnetism, 
however, a new factor enters, namely time. It is also possible that the magnetization 
of the particles with a fairly high coercive force is affected by a fairly weak magnetic 
field, if this field is active during sufficiently long time. This effect becomes 
noticeable during the long geological times and has therefore to be taken into 
account in palaeomagnetic research. Neéel (1949) introduced the conception 
‘relaxation time’ for this phenomenon, which designates the same value as the con- 
ception coercive force. Of a certain ferro-magnetic mineral we may expect that 
a particle with a high coercive force has a high relaxation time too, and reversely. 
After the final formation of the rock up to the moment of collecting, the rock 
specimen has been in the changing Earth field. The direction of the magnetization 
of the particles with short relaxation time easily follows this Earth field and these 
particles cause the disturbing magnetizations in the rocks. But the magnetization 
of the particles with a relatively long relaxation time in respect of the age of the rock 
is only little, or not at all, affected by the present Earth field and hence represents 
the original fossil magnetization. The magnetizations of the particles with the 
relatively short relaxation time are probably characterized by a low coercive force 
and are thus unstable with respect to demagnetization, whereas the original fossil 
magnetization has a high coercive force and cannot easily be removed by demag- 
netization forces. With these assumptions it must be possible to remove the dis- 
turbing magnetizations by partial demagnetization. 


3. Method 


Demagnetization can be performed by an alternating field or by heating, and of 
course by a combination of both. 

In order to demagnetize a specimen by means of alternating fields the sample is 
put in a coil through which an alternating current is sent from the 50c/s mains. 
This coil is placed in a field-free space produced by a Helmholz coil system. 
The field-free space is regularly checked by a milli-oerstedmeter. The maximum 
of the produced alternating field was 300 oersted. For the heating of the 
specimens we had at our disposal a furnace with adjustable temperatures up to 
300 °C. To make the rough specimens suitable for the measurements they were 
cast in orientated position in plaster of Paris in cubes with edges of tocm. This 
system has many advantages; one in particular is that, when the cubes with 
specimens are heated, the plaster acts as an insulating cover and in this way it is 
possible to demagnetize them with alternating fields and to measure a rock 
sample in a hot state. 

The magnetic measurements are executed with the aid of an astatic magneto- 
meter developed by one of the authors. This astatic magnetometer contains an 
astatic system consisting of two small magnets at 5cm mutual distance, suspended 
by a very thin nylon thread. The free period of the system is 20s, the factor of 
astatism is one thousandth, and the sensitivity is 10-* e.m.u. for 1mm deviation 
of the light spot. In the measuring system chosen the magnetic vector of a speci- 
men is measured in three components perpendicular to each other, hereafter 
called A, B and C component. A is the vertical ‘up’ component, B is the south 
component and C is the east component. In order to demagnetize the specimen in 
the alternating current coil, its three components A, B and C are placed one after 
the other parallel to the axis of the coil and so each component is exposed to an 
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increasing and decreasing alternating field in cold or hot state as desired. Immedi- 
ately after this procedure the sample is measured. 'To observe the behaviour of the 
magnetization of the rocks under the process of demagnetization the specimens 
were progressively demagnetized. Thus a sample was treated with ever-increasing 
alternating fields up to 300 oersted, in steps of about 30 oersted, and measured 
after each treatment. Special attention was paid to the alternating current used 
in the coil system. During the experiments of demagnetization by powerful alternat- 
ing fields a disturbing magnetization appeared to be induced in the samples. 
This disturbance, which may be considerable, turned out to be an effect of 
asymmetric distortions in the alternating current. Despite the small distortions 
as caused by higher harmonics (less than 1 per cent) a disturbing field of the 
order of one oersted is produced, when the demagnetization field is 300 oersted. 
The odd harmonics cause a symmetric distortion with regard to the 50 cycles 
current and do not cause any trouble. The even harmonics have to be filtered. 


This was accomplished by a pass filter for the 50 cycles current combined with 
a filter to eliminate the 100 cycles. 


4- Results of the experiments 


The specimens of rocks used for these experiments came from the Estérel, a 
Permian vulcanic region near the coast of the Mediterranean in the south of France. 
This region consists of some dolerite flows, two thin (8m and 15m) rhyolitic 
ignimbrites, one very thick rhyolitic ignimbrite (till 200m) and one rhyolitic 
flow (about 100m, locally called pyromeride). All these flows and ignimbrites are 
often separated by more than 100m thick formations of tuffs, conglomerates, 
arkoses and pelites (Bordet 1951). Of these formations more than 500 oriented 
specimens were collected. A part of these specimens was used in magnetic clean- 
ing tests. Some results of progressive demagnetization are shown in Figures 
1-5. In these examples no corrections for geological strike and dip have been 
made. 

According to the assumption made about magnetic cleaning, it may be expected 
that during the first part of the progressive demagnetization the magnetization with 
low coercive force diminishes and even vanishes without greatly affecting the origi- 
nal fossil magnetization with high coercive force. The consequence is that the 
resulting magnetization changes not only in intensity but also in direction. This 
rotation of the resulting magnetization continues until all the disturbing magnetiza- 
tions are destroyed and only the very stable part remains; it is supposed that this 
is the fossil magnetization. If, after that, more powerful demagnetizing fields 
are applied, only this last-mentioned direction diminishes, which means that the 
direction of the measured magnetization no longer rotates. When the direction of 
magnetization stops moving, even with increasing demagnetizing fields, we may 
assume that the specimen has been cleaned magnetically. This assumption is 
wholly confirmed by the experiments. In some cases the applied field was not 
strong enough to show clearly the last stage of the process of cleaning. However, 
the results obtained up to now may be considered of importance for palaeo- 
magnetic research. 

The most essential features may be seen in Figure 1. This is an example of a 
progressive demagnetization of a dolerite specimen in cold state. The graph shows 
the changes of the three components A, B and C and the total intensity as a func- 
tion of the value of the alternating demagnetization field. The stereographic 
projection shows the calculated directions of the magnetizations. In this and all 
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the other specimens treated by progressive demagnetization, the magnetization 
vector describes a path in the plane defined by the present local Earth field vector 
and the direction of the local permian field. It is remarkable that the B component 
first increases and then decreases with a maximum at about 150 oersted. The 
change at 150 oersted can also be found in the behaviour of the other components, 
This phenomenon clearly points to two magnetizations: one in the Permian direc- 
tion and one in the direction of the present field. It is easy to determine the direc- 
tion of the magnetization that decreases first. The triangle in the stereographic 
projection indicates this direction, agreeing approximately with the direction of 
the present Earth field in the Estérel, which is indicated by an asterisk. The 
bend in the graph of the B- and the other components at 150 oersted is a conse- 
quence of the fact that at this point the original Permian magnetization begins 
to decrease. At this point the disturbing magnetization is not quite removed, 
for next the direction composed of the three components is still moving. Even 
after the sample had been treated in a field of 295 oersted it was not completely 
cleaned. 

The next example is also a dolerite specimen. This sample received the same 
treatment, now at a temperature of 100°C. It seems that with a treatment of 295 
oersted at a temperature of 100°C the magnetic cleaning is finished. This was 
confirmed by treating the specimen with 295 oersted at 150 °C, when only shrink- 
ing of the magnetization vector was observed. Also in this example an agreement 
was found between the direction of the first removed magnetization and the 
direction of the present Earth field. 

Figure 3 shows the progressive demagnetization of a rhyolitic ignimbrite at 
100°C and there it is obvious that the specimen is completely cleaned by a field 
of 295 oersted at 100°C. 

The behaviour of a pelite treated at room temperature is shown in Figure 4. 
Sediments react in the same way as volcanic rocks. The rotation of the direction 
of the magnetization is quite clear. The graph shows that above 150 oersted the 
Permian magnetization is affected but that it decreases slower than does the present 
field magnetization. Again it can be seen that at room temperature a treatment by 
a 295 oersted alternating field is not always enough to complete a magnetic cleaning. 

Finally, Figure 5 presents the magnetic cleaning process of an arkose. In this 
case, fortuitous at 150 oersted, the sample is wholly cleaned. 

The total effect of cleaning a number of samples by a field of 300 oersted at a 
temperature of 150°C is shown in Figure 6. In this figure, data of dolerites, 
rhyolites, pelites and arkoses are collected. At the top left-hand corner the samples 
shown have neither been corrected for geological dip nor by magnetic cleaning. 
At the top right-hand corner the samples have only been magnetically cleaned, 
at the bottom right-hand corner the samples have been corrected both for geological 
dip and by magnetic cleaning. When a mean is taken from the uncleaned and 
cleaned samples there is a difference in direction of about 30° which corresponds 
to a pole displacement of 18°. The individual directions of magnetization of some 
specimens showed rotations up to more than 60° during the process of cleaning. 
(From preliminary experiments on the rhombporphyries of the Oslo region it 
became evident that these deviations may even be greater.*) The magnetic pole 
calculated from 14 demagnetized samples of the Estérel is situated at Long. 144° E 


* Personal communication of Mr R. O. van Everdingen of the Mineralogisch-Geologisch 
Institute, Rijksuniversiteit, Utrecht, who is preparing a doctoral thesis on the palaeomagnetic 
analysis of the rhombporphyries of the Oslo region. 
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Fic. 6.—The effect of cleaning 14 samples of dolerites, rhyolites, pelites and arkoses, 
combined with correction for tectonical strike and dip. 


+ = upper hemisphere 


a.—The measurements before magnetic cleaning and without tectonic corrections. 


O = lower hemisphere 


b.—The directions after cleaning the samples by a field of 300 oersted at a temperature 


of 150 °C. 


c.—The effect of correction for tectonical strike and dip only. 


d.—The directions after magnetic cleaning and after correction for strike and dip. 
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and Lat. 47° N. Of course this position must not be regarded as a definite pole 
position of the Permian rocks of the Estérel, for the research is still in progress 
and the bulk of the specimens has not yet been worked on. 


5. General remarks 


From the well-known formula cot ¢ = }tanJ it can easily be seen that a given 
difference in the inclination of magnetization causes a difference in pole position 
that is a half of the deflection in the inclination of the magnetization when the pole 
distance is about go°, but twice that when the pole distance is about zero. So it 
is clear that even greater errors than have been found can occur in poles derived 
from observed magnetic directions in uncleaned rocks. Even when a great con- 
centration of the directions is obtained, it is not certain that a systematic devia- 
tion is absent. On account of the great concentration, one gets a small cone of 
confidence but cones of confidence do not in fact reveal systematic errors. 

The dipole axis calculated from a direction of magnetization which results 
from two dipole fields is also a resultant of the two corresponding dipole axes. 
This means that the pole position calculated from rocks that are influenced by the 
present field deviates to or from the North. Therefore the derived longitude of the 
pole position will be correct but there may be an error in the latitude of the pole. 


6. Conclusions 


The magnetization of rocks is often the resultant of different magnetizations, 
namely a fossil magnetization that arose during the formation of the rocks, and a 
secondary magnetization which often has the direction of the present Earth field. 
It may be possible that even more than two magnetizations are present but in the 
Permian rocks of the Estérel this has not been observed. 

The disturbing secondary magnetizations are broken down by rather weak 
alternating fields. Generally an alternating field of 300 oersted combined with 
heating up to 150°C is sufficient to remove the disturbing components. The 
original remanent magnetization is clearly affected at room temperature by 
alternating fields stronger than 150 oersted. At higher temperatures this point 
is already reached by weaker alternating fields. However, even after intensive 
demagnetization, this fossil remanent magnetization was well measurable. It 
appeared also to be possible to remove the disturbing magnetizations while saving a 
good deal of the fossil magnetization. 

Finally, from these demagnetization experiments it follows that all pole 
determinations that are based on measurements of remanent magnetizations in 
rocks from which the disturbing magnetization of the present Earth field has not 
been removed should be treated with reserve. 


7. Acknowledgments 


This research has been carried out under the supervision of Professor J. 
Veldkamp, Director of the Geophysical Division of the Royal Netherlands Meteo- 
rological Institute. The authors wish to thank Professor Veldkamp for his valuable 
suggestions concerning the experiments and the interpretation of the measurements, 
as well as for his willingness to read and to criticize the manuscript. 

The authors are also indebted to Professor M. G. Rutten, Director of the 
Geological Institute of the University of Utrecht, who stimulated the palaeo- 
magnetic research and proposed to start with palaeomagnetic investigations in the 
region of the Estérel. 


Bi 


; 
for 
an 
the 
4 
Bye 
C 
I 
3 
I 
] 


Magnetic cleaning of rocks in palaeomagnetic research 319 


Special mention has to be made of Mr R. O. van Everdingen, who co-operated 
in the palaeomagnetic research by his investigations of the rhombporphyries, the 
results of which will be published elsewhere. ; 

Grateful acknowledgment is also due to Mr W. Schuur and Mr G. Klopman 
for their valuable help in preliminary experiments and to Mr E. Ph. F. H. Blokhuis 
and the instrument makers of the Royal Netherlands Meteorological Institute for 
the excellent construction of the instruments. 


Koninklijk Nederlands Meteorologisch Instituut, 
De Bilt, Netherlands. 


1958 July 11 


References 


Bordet, P., 1951. Etude géologique et pétrographique de l’Estérel. Mémoires pour 
servir a l’explication de la carte géologique détaillée de la France, Imprimerie 
Nationale, Paris. 

Clegg, J. A., Almond, M. & Stubbs, P. H. S., 1954. ‘“The remanent magnetism 
of some sedimentary rocks in Britain”, Phil. Mag., 45, 583-598. 

Creer, K. M., 1957a. ‘““The natural remanent magnetization of certain stable 
rocks from Great Britain’, Phil. Trans. Roy. Soe. A, 250, 111-129. 

Creer, K. M., 1957b. ‘“The remanent magnetization of unstable Keuper Marls”, 
Phil. Trans. Roy. Soc. A, 250, 130-143. 

Graham, J. W., 1949. ““The stability and significance of magnetism in sedimentary 
rocks”, J. Geophys. Res., 54, 131-160. 

Graham, J. W., 1955. ‘‘Evidence of Polar Shift since Triassic time”, J. Geophys. 
Res., 60, 329-348. 

Irving, E. & Runcorn, S. K., 1957. ‘‘Analysis of the palaeomagnetism of the © 
Torridonian Sandstone series of N-W Scotland”, Phil. Trans. Roy. Soc. A, 
250, 83-99. 

Nagata, T., 1953. Rock-magnetism, Maruzen Co. Ltd., Tokyo. 

Néel, L. 1949. Annal. Géophys., 5, 99-136. 

Runcorn, S. K., 1952. Handbuch der Physik. Springer Verlag, Berlin, Vol. XLVII, 
PP- 57°-597- 

Rutten, M. G., van Everdingen, R. O. & Zijderveld, J. D. A., 1957. ‘‘Palaeomag- 
netism in the Permian of the Oslo Graben (Norway) and of the Estérel 
(France), Geologie en Mijnbouw, 19, 193-195. 


| 
f 
S 
ie 
L 
t 
5 


A Magnetic Survey off the West Coast of the United States 
between Latitudes 32° and 36° N, Longitudes 121° 
and 128° W 


R. G. Mason 
(Received 1958 July 30) 


Summary 

A total magnetic intensity map covering the area between latitudes 
32° and 36° N and longitudes 121° and 128° W is presented. It shows a 
series of narrow anomalies of about 400 y amplitude trending north— 
south for 300 miles over a substantially flat abyssal plain. Some possible 
interpretations are discussed. The anomalies arise at shallow depths 
beneath the ocean bottom and are consistent with the presence of 
relatively thin layers of magnetized material overlying a less magnetic 
crust; they could equally well arise from topography of the upper surface 
of a continuous magnetic basement for the most part concealed by an 
overlying blanket of sediments. 

The anomalies are intersected by the Murray fracture zone, which 
interrupts their pattern in such a way as to suggest a right lateral dis- 
placement of about 84 nautical miles. 


1. Introduction 


Since the latter part of 1952 the Scripps Institution of Oceanography has been 
towing a total field magnetometer astern one or other of its research ships engaged 
in routine oceanographic investigations. With a few exceptions, as for example 
in the study of the Tonga Trench (Raitt, Fisher & Mason 1955), the results have 
been restricted to long and unrelated profiles obtained during passage of the ship. 
These have limited application since they give no indication of the lateral extent 
or direction of anomalous magnetic trends and therefore provide no basis for 
quantitative geological interpretation. 

In August 1955 the opportunity arose to tow the magnetometer behind the 
U.S. Coast and Geodetic Survey Ship Pioneer engaged in making detailed hydro- 
graphic surveys with closely spaced and accurately surveyed lines. As a result an 
extensive area of the deep Pacific off the west coast of the United States (Figure 1) 
has now been mapped with the detail and accuracy generally associated with 
airborne magnetometer surveys over the land, and advantage can be taken of the 
techniques which have been developed for the interpretation of such surveys. 

This paper deals with the practical aspects of the survey and discusses some 
possible interpretations consistent with present day speculation on the nature of 
the ocean floor. The work was supported by the U.S. Office of Naval Research. 


320 


“atts 


tes 


A magnetic survey off the west coast of the United States 321 


2. The survey 


A notable feature of the field operation was the complete coverage of the area 
with a regular grid of accurately surveyed lines, usually only five miles apart. 
The position of the ship was determined at approximately eight-mile intervals 
by means of a radio navigation system using ground wave propagation and fixed 
beacons on shore. The probable error of an observation varies with time of day 
and distance from the beacons but is of the order one-tenth of a mile. The ship’s 
tracks were for the most part directed east-west, that is across the principal grain 
of the anomalous magnetic field, but sufficient lines were run in the north-south 
direction to eliminate any possibility of ambiguity in contouring. 


SAN 
FRANCISCO 


35°N 


30°N 
130° W 125°W 120°W 115°W 


Fic. 1.—Map showing the area covered by the survey. 


The magnetometer was an ASQ-3A MAD, originally designed for detection 
of submarines, modified for geophysical investigations according to plans developed 
by the U.S. Naval Ordnance Laboratory and further modified to improve its long 
term stability and render it suitable of use over long periods at sea. It is a gimbal 
mounted flux gate oriented in the direction of the Earth’s magnetic field by two 
other flux gate elements which control orienting servo motors. The unit compris- 
ing the gimballed flux gates and associated motors was housed in a streamlined 
fish of plastic material which was connected to the electronic and recording equip- 
ment on the ship by a 500 ft electrical tow-cable. At that distance the effect of 
the ship on the measured field was negligible, ranging between about +5 y, 
depending upon the ship’s heading. 

The magnetometer was calibrated ashore against a proton precession magneto- 
meter at approximately monthly intervals. Its total drift over the period of the 
survey was 20y. ‘Taking into account the principal sources of error, namely 
inherent drift, inaccurate adjustment, and the effect of the ship, the probable 
absolute error of a single observation is about 10 y. With proper adjustment the 
relative error between any two points on the same line is negligible. 
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The magnetic map (Figure 2) was constructed by first plotting the measured 
field value, in 50 y units, along the ship’s track and contouring the result with an 
interval of 50 y. The normal geographic change of the total magnetic intensity 
was determined by dividing the map into 20 mile squares, computing the mean 
field for each square and smoothing these values in both the north-south and 
east-west directions by a least squares process. A plot was then made of these 
smoothed values and this was subtracted from the original plot and the difference 
contoured with an interval of 50 y. 

No attempt was made to remove the usual time variations of the Earth’s 
magnetic field from the data. However, records of the total field at Palo Alto, 
near San Francisco, were available for several periods of a few days around the 
time of the survey and these were found to agree closely with the corresponding 
records from Tucson, Arizona, magnetic observatory. Since Palo Alto is about 
700 miles WNW of Tucson but less than 100 miles to the north of the area under 
investigation it may be concluded that Tucson records give a fair indication of 
time variations over the area of the survey. There were no magnetic storms during 
the course of the survey and the extreme range of the mean hourly values of the 
total field at Tucson was only 60 y for the whole period. The mean daily range of 
the hourly values was 25 y. These variations are small compared with the ampli- 
tudes of the principal anomalies and although they may introduce small errors in 
the pattern of the contours, particularly in the areas of shallow magnetic relief, 


they can have little effect on the positions and magnitudes of the more significant 
anomalies. 


3. Preliminary interpretation 

The magnetic data and the one relevant seismic section at that time available 
were sent to the Gravity Meter Exploration Company of Houston, Texas, for 
preliminary interpretation. On the basis of a number of depth estimates, using 
techniques which they describe as ‘“‘such as might be devised putting into effect 
the principles outlined in GSA Memoir No. 47” (Vacquier & others 1951), they 
reach the general conclusion that the anomalies arise from relatively shallow struc- 
tures. Referring in particular to the observed anomaly across the section marked 
A-A’ in Figure 2, which corresponds approximately to the location of the seismic 
section of Figure 3, their report states, in paraphrase: 


“In an attempt to determine the size and shape of the structures responsible 
for the anomalies a number of theoretical anomaly curves were computed. From 
Figure 2 it can be seen that a number of the more pronounced anomalies, parti- 
cularly in the western part of the survey, are oriented north-south geographic. 
Figure 4 shows the computed anomaly for an east-west crossing of an infinite 
north-south block. The thickness of the block is 0-4 times the depth of burial, 
or approximately 2 km. The susceptibility value used is quite high, 0-o109 c.g.s. 
The asymmetry of the computed anomaly arises from the difference between 
the orientation of the block and the orientation of magnetic north. Superimposed 
on the computed anomaly is the measured anomaly across section A-A’. The 
agreement is striking. We feel that the anomalies are caused by magnetized material 
in the form of a relatively thin slab above the normal oceanic crust. This would 
correspond with the layer marked ‘volcanics’ on the seismic section. 

A point of interest is the high susceptibility value that is necessary for the 
theoretical curve to fit the observed data. It is high even for the basic igneous 
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materials that might be expected in this area. It is our opinion that such a high 
value may be indicative of appreciable remanent magnetization of the slabs 
obtained during cooling in the presence of the Earth’s field.” 


4. Further two-dimensional solutions 


The suggestion that the anomalies arise from lava flows is consistent with 
contemporary speculations on the nature of the ocean floor (Revelle & Raitt 
1957). However, it is not necessarily the only possible explanation. The remark- 
able regularity of the more prominent anomalies throughout the several hundred 
miles of their length covered by the survey points to a relatively simple cause and 
encourages a more thorough investigation of possible two-dimensional solutions. 

Referring again to magnetic profile A~A’, the solution is restricted in vertical 
extent by the floor of the ocean at about 4-6 km and, presumably, by the Mohoro- 
vicic discontinuity which, according to the seismic data, is here at a depth of about 
11 km. It is further restricted in that physical considerations impose a maximum 
value on the depth of its upper boundary, as will be shown later. If certain simpli- 
fying assumptions are made about its shape and state of magnetization, then 
with a suitable choice of depth and susceptibility it is possible to derive a 
unique section, within the above limits of depth, which will fit the anomaly 
exactly. 

Figure 5a shows one such section, derived on the basis of a uniformly magne- 
tized single continuous section with a plane lower surface at a depth of 6-3 km, 
corresponding to the base of the “volcanics”, and with induced and permanent 
magnetization in the direction of the present field equivalent to a susceptibility 
of 0-015. The derived upper surface, which is unique, extends at its highest point 
into the sedimentary layer. This and the sections discussed later were built up, 
by trial and error, of several hundred infinite rectangular strips each 0-5 km wide 
by o-2 km thick. The observed curve was matched to within one or two y. 

Similar computations have been made for a plane upper surface. Figure 5b 
shows three such sections based respectively on equivalent susceptibilities of 
0°015, 0°0075 and 0-005, the upper surface in each case lying at a depth of 5 km, 
corresponding to the base of the sediments. The first is contained almost entirely 
within the ‘“‘volcanics” layer, the second extends well down into the crust and has 
steeper edges, the third is terminated by the Mohorovicié discontinuity. In 
general, for a fixed upper surface an increase in susceptibility leads to a thinner 
section of similar shape but more drawn out at the edges, a reduction makes the 
section thicker and its sides steeper. In theory there is no upper limit to the suscep- 
tibility, and corresponding thinness, of the section, but the depth limitation im- 
posed by the Mohorovicic discontinuity leads to a lower limit of 0-005. 

The result of increasing the depth of the upper surface is illustrated in Figure 
5c. The sides become increasingly steeper until at a depth of 6-2 km, for a suscep- 
tibility of o-o15, the structure becomes an almost rectangular block with its lower 
surface at 7°38 km. This block is the deepest structure of susceptibility 0-015 
that can be fitted to the anomaly, though it is possible to make the upper surface 
still deeper, and the lower surface correspondingly shallower, by increasing the 
susceptibility. The limiting case is an infinitely thin sheet, of infinite suscepti- 
bility, at a depth of 7-0 km. 

In their report, the Gravity Meter Exploration Company draw attention to 
the fact that several of the north-south anomalies appear to have two maxima, 
with an intervening minimum, and suggest that the lines of minimum field have 
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some structural significance. Figure 6 shows the magnetic profile along the line 
marked B-B’ in Figure 2, which crosses one such anomaly, together with a derived 
section based on the assumptions of a plane upper surface at a depth of 5 km and 
uniform magnetization equivalent to a susceptibility of o-o15 in the present direc- 
tion of the Earth’s magnetic field. In general this section responds to changes in 
depth and susceptibility in the same way as the sections of Figure 5, except that 
for greater depths the section breaks into two isolated parts. On the assumption 
of a depth of 11 km to the Mohorovici¢ discontinuity, the least possible suscepti- 
bility is o-or. 
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Fic. 6.—Two-dimensional interpretation of magnetic profile B—B’ for a 
field of 50000 y, dip 59°, declination 17°E, susceptibility 0-o15. 


5. Remanent magnetization 


Attention was drawn earlier to the high susceptibility values needed to explain 
the observed anomalies. The preceding analysis shows that for section A-—A’ 
the least possible susceptibility is 0-005, corresponding to the improbable case of 
a structure extending from the ocean bottom to the base of the crust. For the more 
plausible structures the susceptibility must be well in excess of 0-o1, or higher 
than the known susceptibilities of any of the common rocks. This strongly suggests 
that the rocks involved here are basic igneous rocks, which are characterized not 
only by a high susceptibility but also by a high intensity of remanent magnetiza- 
tion, commonly several times greater than the intensity of induced magnetization 
in the Earth’s field. Thus in the present case, where the Earth’s field is 50 000 y, 
an apparent susceptibility of 0-015 could be accounted for by a true susceptibility 
of 0-005 together with an intensity of remanent magnetization of 0-005 G in the 
direction of the Earth’s field, both acceptable values for a basic igneous rock. 
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The interesting question arises as to whether, since the direction of remanent 
magnetization will correspond to the direction of the Earth’s field at the time when 
the rocks were formed, in all probability different from the present direction, 
anything can be learned about the direction of magnetization that may throw 
light on the age of the formations. Two factors bear on the problem. Geological 
evidence (Menard 1955) points to a Cretaceous or early Tertiary age for the 
Murray fracture, and the magnetic structures are older than this since their pat- 
tern is interrupted by it. Palaeomagnetic measurements (see Runcorn 1956 for 
a summary of relevant information) suggest that since Jurassic times the direction 
of the Earth’s magnetic field in North America has centred around that of a dipole 
field, with normal or reversed polarity, but that during earlier periods it was 
quite different. 

Figure 7 shows four sections computed to test the effect of varying the horizontal 
direction of remanent magnetization. The susceptibility was taken to be 0-005 
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Fic. 7.—Two-dimensional interpretation of magnetic profile B—B’ for a 

field of 50000 y, dip 59°, declination 17° E, susceptibility o-oo5 and inten- 

sity of remanent magnetization 0-005 G dipping at 59° in azimuths 
(a) 34° W, (6) 0°, (c) 17°E, (d) 34°E. 


in all cases and the intensity of remanent magnetization 0-005 G dipping at an 
angle of 59°, the present dip, but in four different azimuths, namely 34° W, 0°, 
17°E (the present magnetic declination) and 34°E respectively. Although the 
section corresponding to the present direction is perhaps the most symmetrical, 
there is no compelling reason for preferring one to the other or for rejecting any 
one. This lack of sensitivity to azimuth is a consequence of the fact that only the 
components of magnetization perpendicular to the strike, i.e. the vertical and the 
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east-west, contribute to the anomaly, and for steep dips the former predominates. 
For a dip of 70° all computed sections would lie between the limiting sections of 
Figure 6, irrespective of azimuth. A more thorough investigation is beyond the 
scope of this paper. However, it seems clear that, since any remanent magnetiza- 
tion is required to augment rather than to oppose the present field, the vertical 
component must be substantially downwards and this makes any appreciable 
degree of reversed magnetization unlikely if the magnetic poles were anywhere 
near their present positions. 


6. Conclusions 


Although the discussion of the preceding sections was based upon a much 
simplified picture of the situation, there can be little doubt that the major positive 
anomalies are underlain by an excess, probably totalling at least 1 km in thickness, 
of a material highly magnetic by comparison with adjacent formations and almost 
certainly a basic igneous rock. Its upper surface is relatively shallow and probably 
lies in the ‘“‘volcanics” layer of the seismic section, or at most no deeper than in 
the upper part of the crust. Its extent in depth is uncertain but it may be supposed 
to terminate at a depth no greater than the Mohorovicic discontinuity. Despite 
its thickness it shows little or no correlation with topography, as far as can be 
judged from the best available topographic maps, suggesting that it is either 
buried in consolidated sediments of relatively high seismic velocity or emplaced 
in less magnetic igneous rocks. 

The most obvious and acceptable geological explanation is that the excess 
material represents lava flows which have spread out over the floor of the ocean, 
possibly filling pre-existing troughs, or depressions produced by the weight of 
the accumulating flows. This picture is consistent with the section of Figure 4a 
and the upper section of Figure 4b, which thin out at the edges in a manner charac- 
teristic of lava flows. Within the limits of the survey the major anomalies cover 
an area of about 50000 sq km. In extent therefore, and in thickness, the flows 
would rank with the largest continental lava fields, such as the Columbia River 
basalts which cover an area of 120000 sq km and in places reach thicknesses 
approaching 2 km. An alternative hypothesis, consistent with the section of Figure 
4a, is that the excess material represents topography of the upper surface of the 
crust, itself supposedly the more magnetic layer, brought about by marginal 
faulting or by folding. In that case it must be presumed that the topography has 
for the most part been covered by a blanket of sediments. 

Whatever their origin and form, the rock bodies which give rise to the north— 
south anomalies indicate lines of weakness in the ocean floor which were active 
prior to the formation of the Murray fracture zone. An offset of the pattern of the 
magnetic features as they cross it suggest a right lateral displacement of 84 nautical 
miles along the Murray fracture. This displacement, the north-south lineation of 
the anomalies and the right lateral displacement along the San Andreas fault, 
may well be related by a single structural hypothesis, possibly connected with a 
westerly movement of the North American continent relative to the floor of the 
Pacific. 

This paper has been mainly concerned with the north-south anomalies because 
they are the most striking feature of the magnetic map and because the seismic 
section provides a framework for discussion which is lacking in the case of the 

NNE-SSW trending anomalies on the eastern portion of the map. It is hoped 
that further geophysical investigations, particularly seismic, planned on the basis 
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of this map will help to throw more light on the complicated structural situation 


in this area and hence on the more general problem of the relationship between 
continent and ocean. 
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Summary 


Gravity observations made with gravity meters over large gravity 
ranges have shown the instruments to be less satisfactory than is 
desirable. Figures are given showing the errors which have been 
found using three gravity meters under various field procedures and 
some possible causes of error are discussed. The scale constants of 
the gravity meters are shown to change substantially with time, and a 
system is described which was built into one gravity meter in order to 
check the scale constant at any time. The results obtained from this 
device are tabulated and discussed; no satisfactory mechanism has 
been produced to explain the change in scale-constant, and the 
calibration device has some use, although not so much as had been 
hoped. 

It is concluded that for the accurate measurement of large gravity 
differences careful operation of the gravity meters is essential, the 
calibration of the instrument being checked regularly; and that pendu- 
lum measurements may be preferable for the measurement of very large 
gravity differences. 


1. Introduction 


The Dominion Observatory makes regular long-range gravity meter observa- 
tions. Wherever possible, connections are made to any available pendulum 
bases, in order to improve the accuracy of the National Gravity base network. 
The connections to pendulum bases also afford a means of calibrating the instru- 
ments, the accuracy of calibration depending on the reliability of the pendulum 
values and the gravity difference between the bases. 

Evidence has been accumulating that the scale constants of the gravity meters 
change with time (Innes & Thompson 1953). Accordingly two convenient pen- 
dulum sites, at Ottawa and Washington, were chosen to terminate a calibration 
line, over which the gravity meters could be checked regularly (Innes 1958). 
More detailed information on the behaviour of the gravity meters which has been 
collected over the past three years, using this calibration line, has prompted the 
construction of a device for quickly checking the calibration of one of the North 
American instruments. The purpose of this paper is to describe both the behaviour 
of the instruments and the calibration device. 


* Published with the permission of the Deputy Minister, Department of Mines and Technical 
Surveys, Ottawa, Canada. 
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2. Methods of gravity measurement 


The Dominion Observatory makes an invariable practice of following a looping 
procedure to eliminate drift when calibrating gravity meters and when establishing 
primary control points. Drift curves are then plotted so that errors due to drift 
are minimized. Aircraft are used for long-distance surveys and cars for shorter 
distances. The distance between stations is arranged so that a reading is taken 
every hour or less. The length of a “‘leg” is greater when an aircraft is used, but 
it is found that the accuracy of the reading for one “‘leg”’ is less than when a car 
is used. However, over a given distance, the increased number of “legs” involved 
reduces the accuracy of the car work, so that it is comparable to that done by air. 

The measurements between Ottawa and Washington have been made using a 
car exclusively. Three gravity meters have been used—a Geodetic type North 
American instrument (No. 137), a normal North American instrument (No. 85), 
and a Worden gravity meter (W. 44). These instruments will generally be identi- 
fied by number henceforth. The distance between the two bases is about 700 
miles, and twenty-four intermediate stations are used. The route has changed 
slightly since the programme was started, as have some of the stations. However, 
the changes have been made so that never more than two consecutive stations 
were altered, and in all cases the old and new stations have been connected by 
looping. 


3- Results of base-line looping 


A tabulation of all the readings made during the many runs would be tedious, 
and the important facts may be presented more simply. 

The accuracy of repetition depends to a considerable extent upon the instru- 
ment involved. The calibration line stations have been measured so many times 
that, assuming values of gravity at Ottawa and Washington, the probable error of 
the adopted gravity difference between adjacent stations is about 0o-o1mgal. If 
the scale constants of the instruments are adjusted to give the best fit, for each 
traverse, then W. 44 has an r.m.s. error of 0-04mgal, No. 137 one of 0-05 mgal, 
and No. 85 has an r.m.s. error of 0-o8mgal for the gravity difference between ad- 
jacent stations. It should be noted that the Worden instrument, W. 44, is used 
with the smal] dial only; this instrument has a Geodetic dial, and if it is used the 
error is approximately tripled. The scale constants of the instruments are cal- 
culated to have accuracies of 1 part in 10000 or better for a given adopted gravity 
difference between Ottawa and Washington. The variations of scale constants 
of the three instruments are plotted against time in Figures 1 to 3, and it is evident 
that there are real changes in the calibration. It has been suggested that Worden 
instruments have a constant whose value changes with temperature by approxi- 
mately 1 part in 1000 for 17°C, corresponding to the temperature coefficient ot 
the shear modulus of elasticity for fused silica. Variations of scale constant deter- 
mined with W. 44 cannot be explained in such a manner. 

It can be seen that the scale constants of all the instruments change when they 
are repaired, as might be expected. The thermostat relays in the North American 
instruments occasionally weld themselves in a closed position, and the temperature 
of the instrument rises exponentially towards about 100°C. When the instrument 
has been overheated in this way, the scale constant usually changes, in one case by 
more than 1 per cent. There have been no facilities for measuring either the dura- 
tion of the incident or the maximum temperature reached, although it is known that 
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Fic. 1.—Worden No. 44: variation of scale constant with time. 
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the temperature has usually considerably exceeded 50°C. As well as the above 
two changes, of which the cause is known and can be guarded against, other changes 
occur for no known reason. The Worden instrument in particular, which is 
constructed entirely of fused silica, would be expected to have a fixed calibration, 
but that this is not so can be seen from Figure 1. Whilst calibrating other gravity 
meters recently it was found that the scale constants differed depending on whether 
the calibration was made in a direction of increasing or decreasing gravity. In 
the case of three Worden instruments, the variation ranged from 1 part in 1 000 
to an amount which was less than the error of determination. This effect has yet 
to be investigated more fully. 


4. Causes of instrumental error 


The Worden gravity meter, which has no thermostat, shows a very strong 
correlation between temperature change and drift rate. If W. 44 is mounted in a 
thermostat, the drift curves are very much improved. When the drift curves of 
the two North American instruments are compared, they show many similarities. 
In order to check that the thermostats were working correctly on the two North 
American gravity meters, thermometers were mounted on the case to read the 
temperature of the lid of the inner temperature enclosure to o-o1 deg C. These 
have shown a variation in the inner temperature of up to 0-05 deg C during the 
day, but there appears to be but little connection between the temperature and the 
drift curve. The barometric pressure at the stations does not appear to influence 
the drift directly, but changes of barometric pressure may affect the temperature 
of the gravity meter, the effect being integrated over a considerable period. The 
frame in which the gravity meters are carried in the car is padded and only extends 
about half-way up the meter, whilst the thermometer mentioned above is mounted 
close to the thermostat element at the top of the instrument; it is possible that a 
thermal gradient may be set up in the meter, its value depending on the external 
temperature. This effect would also be integrated. It is noticeable that on a day 
where the temperature stays fairly steady, the drift tends to be small. In this 
connection it may be mentioned that whilst Martin publishes results which point 
out the desirability of correcting for the internal temperature (Martin 1954, pp. 
29-30), his initial readings at a station, which correspond to normal looping read- 
ings, are not improved by the corrections. His method of reading the instrument 
over a period of around two hours is not suitable for normal survey work in Canada, 
as many gravity stations are far from suitable buildings and the temperature inside 
a car often changes very rapidly when the car is stopped for some time. 

The gravity meters were taken from Ottawa to Washington in May 1956, 
reading at the intermediate stations without looping. The journey occupied two 
days and one night; the over-night drift was allowed for, and then a uniform 
drift rate postulated over the two days. Tidal effects were also computed. The 
values of gravity at the intermediate stations are known with an accuracy of 0-o1 
mgal from previous looping, so the error at each station can be found. No. 137 
gave the best results, with a r.m.s. error of 0-13 mgal at the intermediate stations, 
whilst No. 85 and W. 44 gave errors of 0-40 and 0-52 mgal respectively. Compared 
with the looping procedure carried out later that month, the error for No. 137 is 
doubled, that for No. 85 is quadrupled and that for W. 44 is multiplied fourteen 
times. The calibration constants from the run differed from that obtained by loop- 
ing by 1 part in 1000 to 1 part in 200. The results from the two North American 
instruments show how the observations taken without Jooping emphasize the defect 
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of a slightly inferior instrument, whilst the particular Worden used is seen to be 
unsuitable for this type of work. 


5. The calibration device 


When long-range measurements are made, the changes in calibration of the 
instruments are of importance. An attempt was made to build a device to check 
the variation of the scale constant of No. 137, by loading the beam of the gravity 
meter with a weight and measuring the change in reading. A small sapphire ball 
is used as the weight and is captive in a metal cage. A carrier is mounted on the 
beam of the gravity meter and has four arms which project vertically. When 
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Fic. 4.—Cage assembly. Fic. 5.—Ball carrier. 


the cage is lowered, the ball rests on these arms, whilst the cage does not touch 
either the carrier or the ball. Figure 4 shows the cage, ball and carrier in the posi- 
tion where the ball is resting on the carrier and Figure 5 shows the construction 
of the carrier, which is mounted on the gravity meter beam. Consecutive readings 
with the ball loaded on and removed from the beam agree within o-1 division, 


which is as well as the instrument can be usefully read. The change in reading is 
about 350 divisions. 


6. The effect of changes of gravity and changes of mass on the North 
American gravity meter 


The working system of the North American gravity meter is shown in Figure 6. 


The reset control of the gravity meter moves the top end of the reset spring 


approximately vertically, whilst the reading system moves the fulcrum ‘“‘O” 
vertically. 
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Taking moments about ‘‘O” when the value of g is go, and the beam is hori- 
zontal, 


Mgoc sin@ = Td+7b sind, 


M being the total mass of the beam, c the distance from the fulcrum to the centre 
of gravity, T the tension in the main spring and 7 the tension in the reset spring. 


Reset Spring 


Reset Spring Maine, Reset 
Spring Spring 


(8a) M (g, +6 g) 
IM (g,+8g) 


Fic. 6.—Principle of the North American gravity meter. 


The main spring has “proportional length” properties—that is 


T = kl. 
Now 
l b 
— = —, and d=asinB; 
sinB 
therefore ‘ 
bsind 
Td=k asin B = kab sin6; 
sinB 
therefore 


Mgoc sin = kab sin @+7b sind. (1) 


If the instrument is not reset, and the main spring length is kept constant, then when 
the value of g is g9 +g, the fulcrum must be lowered by an amount da to maintain 
equilibrium. 
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The equation then becomes 


M(go+ dg)c sin 6” = k(a+8a)b sin + 7b sing”. 


Now 
da 
b+ 86 = 
therefore 
sing’ = sin¢ cos cos 
and 


sin ¢ cos(Sa/b) cos ¢ sin(Sa/b) 
: + ‘ 
sin sin 0” 


8a\2 
Now @=90, and siné” =1 to the second order. 


Similarly, 
5a being small compared with b. 


8 
Therefore  Mego+Mcdg = kba+kb3a+7b sind-+7b cos 


Subtracting (1), 


and 


da Mc 


It may thus be seen that, unless ¢ is go°, the scale sensitivity is dependent on + 
to some extent. An estimate of the magnitude of this effect can be made. 

The value of g, when Sa = 0 and + = 0 is of the order of 975 gal. When gy 
= 976 gal, 


kb+r7cos¢ 


7b sing = Mgoc — kab (from equation 1). 
At &=81, Mgic = kab (3) 
therefore 
tb sing = Mc(go—81) 
and 
(4) 
Mc btan¢g 


In the North American instrument, the various quantities have values of the 
order quoted below: 


grammes 65—10 cm 


c—10 cm ¢—75° 
cm. 
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Mcig = kbia+7 cos pda 
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(3) 


(4) 
the 


From (3), 


From (4), 


— 
Mc btandé 10x4 


Thus the scale constant will change by an amount of the order of 1 part in 4.000 
for a change in g of 1 gal. It should be noted that this is so only when the instru- 
ment is reset. 

When the ball is lowered on the beam and the instrument readjusted, initial 
conditions being as given in equation (1), 


(M+68M)goc sin 6” = k(a+8a)b sin + 7b sing”. 
Combining this with (1) as before, we get 
5M kb+7cos¢ 
Multiplying this by g9/M 
kb+7 cos¢ 
which is equal to 5g/5a, the calibration constant. Thus, in order to convert the ball 


“reading” to scale constant, correction should be made for the value of g where the 
test is performed. 


7. Results of the calibration device 


The instrument is more sensitive to levelling when the ball is loaded on the 
beam—presumably because the centre of gravity of the beam has been raised 
slightly. This might cause an error of possibly 1 part in 2 000, if the instrument is 
not correctly levelled. but it should be well below this figure normally. A periodic 
error in the screw or dial may adversely affect the accuracy—it is thought that the 
error due to this cause may be about 1 part in 3000. The periodic error will of 
course affect all measurements, to an amount of up to 0-03 mgal. 

The calibrating device has shown the scale of the instrument to be somewhat 
non-linear, as is shown in Figure 7. This amount of variation is to be expected 
from a fairly high-quality micrometer screw. The scatter between the sets of 
observations made during a single month is probably due mainly to periodic 
errors in each turn of the screw and dial. Comparisons of the complete runs with 
the ball are made by plotting curves as in Figure 7 and estimating the change in 
constant, and also by comparing the mean values of the difference in reading caused 
by the ball on each occasion. The two methods always agree within 0-1 division. 

Table 1 gives the mean values of the several ball tests performed on the instru- 
ment, the values of g at the places where the tests were performed, the value of the 
ball tests divided by g (which should be proportional to the scale constant), 
occasions on which the instrument overheated, and occasions when the scale con- 
stant was determined along the Ottawa-Washington line. Investigations into the 
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effect of levelling the instrument showed that the ball test gave results 0-15 divi- 
sion different depending on whether the instrument was levelled when reading 
5 000 divisions or o divisions; accordingly, tests made since June 1957 were per- 
formed with the instrument levelled at around 2 500 divisions. The instrument is 
level when the axis joining the fulcrum “‘O”’ and the centre of gravity of the mass 
is horizontal; as the fulcrum moves with a change in reading, so the “level” posi- 
tion changes slightly. 


348-0(- 


~. September 1956 


3 46-0 ! J 
1,000 2,000 3 4,000 5 ,;000 


Scale Reading, Divisions 


Fic. 7.—Results from Calibration Device. 
Ordinate: Effect of Ball, scale division. 


+---------- + Washington 
© © Ottawa 
e—--—--—-+-—e Syracuse 


It is regrettable that, during the past two years, the instrument overheated 
on an unusually large number of occasions. There is a tendency for the first ball 
test performed after recent overheating to give a high value in column 5, and for 
the rest of the values to agree within much better than one part in 1000. The 
value obtained at Resolute is exceptional in that the variation with scale-reading 
was quite different from normal, showing very small variations with scale reading. 
From the internal consistency of reading, and accuracy of repetition, the accuracy 
of the test would seem to be of the order of one part in 3 000. 

When the values given by the ball test are compared with those from calibration 
runs, a very different picture may be seen. The change in constant from May 1956 
to September 1956 was 0-29 per cent, whilst the ball test showed a change of o-11 
per cent. Between September 1956, and June 1957, the scale constant changed 
by 0°55 per cent, whilst the ball test indicated a change of o-go per cent. From 
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Table 
Results of ball test on North American Meter No. 137 


Effect of Value of Effect of 


Date Place ball gravity ball + g Comments 
1956 div. gal. 
May Ottawa 346°73 980-62 035358 
May Syracuse 346-6 980-40 0°35353 
May Calibration constant 
0°24116 
July Overheated for short 
period 
September Ottawa 347°15 980-62 0°35402 
September Syracuse 347°24 980-40 0°35418 
September Washington 347°2 980°12 0°35424 
September Calibration constant 
0°24045 
1957 
January Overheated for con- 
siderable time. 
February Binghamton 350°28 980-26 0°35733 
April Resolute 351°65 982-86 0°35788 
April Churchill 350°97 981°77 0°35748 
April Ottawa 350°5 980-62 035742 
June Overheated for short 
period 
June 10 Ottawa 350°7 980-62 035763 
June 10 Syracuse 350°35 980-40 0°35735 
June 10 Calibration constant 
0°23914 
June 23 Manitoulin 350°3 98065 0°35721 
June 23 Manitoulin Overheated 
June 24 Manitoulin 350°58 980:65 0°35750 
July Swift Current 350°17 980-88 0°35700 
July Calgary 350°19 980°83 0°35703 
August Prince Albert 350°4 981-23 0°35710 
August 28 Overheated 
August 29 Saskatoon 350°8 981-14 0°35754 
September Ottawa 350°17 980-62 0°35709 
September Syracuse 350°06 980-44 0°35704 
September Calibration constant 
0°23901 


June 1957 to September 1957, the constant decreased by 0-05 per cent, whilst 
the ball test indicated no significant change. The discrepancies here are far greater 
than expected from experimental errors and may possibly give an indication of the 
mechanism by which the scale constant changes. The ball test may be useful for 
checking small changes in scale constant and is certainly useful for checking the 
time at which the changes occur. For large changes in the constant, the accuracy 
is so low as to be virtually useless. 

No satisfactory mechanism has been produced to explain the changes in scale 
constant; the fact that the scale constant can change by 1 per cent with a change in 
the reading of much less than o-1 per cent indicates that the cause is probably not in 
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the main spring. The disagreement between the ball test and calibration would 
seem to indicate that the mathematical analysis given above is too much simplified, 
To produce an error of 0-4 per cent, the position of the ball would have to be moved 
along the beam by about 0-020 in, which is unlikely, nor is the mass of the ball 
likely to change by 0-4 per cent. 


8. Conclusions 


The calibrating device appears to have some use, although a very limited one. 
The behaviour of Geodetic instruments leaves much to be desired, especially 
with regard to temperature changes—which are apt to be greater with air work than 
with car work. For gravity measurements of the highest accuracy, a looping 
procedure should be carried out—the errors previously attributed to ‘‘tares” may 
well have a thermal origin. The calibration constant of Geodetic meters needs 
checking at least once a year, and preferably every three months. Due to the 
uncertainty in the calibration constant, it would seem likely that good pendulum 
observations may be preferable to gravity meter observations when gravity changes 
of the order of 1 gal are being measured; in addition, the behaviour of instruments 


used to interpolate between such pendulum stations needs to be thoroughly in- 
vestigated. 


Dominion Observatory, 
Ottawa, 
Ontario, Canada: 
1958 August 5 
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Research Note 


Determination of the Earth’s Gravitational Potential 
from Observations on Sputnik 2 (19578)* 


A. H. Cook 


(Received 1958 September 24) 


Summary 


R.A.E. observations of the rate of rotation of the orbital plane of 
Sputnik 2 are combined with Jeffreys’s (1943, 1952) equations for the 
coefficients J and D in the spherical harmonic expansion of the Earth’s 
gravitational field to give revised estimates of J and D. The relative 
uncertainties of the results emphasize the great value of satellite 
observations and the need for a revision of the spherical harmonic 
analysis of surface gravity. Comparison of the value of e~! with that 
deduced from the constant of precession indicates that the Earth is not 
in a hydrostatic state. 


Merson & King-Hele (1958) have shown that kinetheodolite observations of 
Sputnik 2 give a rate of rotation of the orbital plane significantly different from 
that calculated from Jeffreys’s (1952) values for the coefficients in the expansion 
of the Earth’s gravitational potential. The precision of this observation exceeds 
that of Jeffreys’s values but the coefficients of the second and fourth harmonics 
cannot be determined separately with the existing data. 

Let the gravitational potential be written as usual in the form 


(=-2 IPA cos 8) + Paco a). 

Here f is the Newtonian constant, M the mass of the Earth, R the Earth’s 
equatorial radius 6 378-10km and (r, 6, ¢) spherical polar coordinates. All other 
terms in the expansion are ignored since the only other one which the satellite 
observations could reveal is that proportional to P3 (cos 6) and the evidence indi- 
cates that this is very small. 

King-Hele (1958) has shown that with such a potential, the orbital plane of a 
satellite rotates at a speed 
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Here g, is the equatorial value of gravity, 7 is the harmonic mean radius and a 
is the semi-major axis of the satellite orbit, and « is the inclination of the plane of 
orbit to the equator. According to King-Hele’s analysis, the terms depending on 
J and D in the motion of the satellite are separable from those depending on air 
resistance. Should an interaction between these terms have been omitted there 
would be a systematic error in this result but, since substantially the same value 
of the ratio of the observed value of Q to that calculated from Jeffreys’s values of 
J and D was obtained at different periods in the life of the orbit, it seems that 
such an error is unimportant. 

Let SJ and 8D be differences from Jeffreys’s values namely, J = 1-637 x 107% 
and D = 10-6 x10-$, and let SQ. be the difference between the observed value of 
. and that calculated from Jeffreys’s values of J and D. 

Writing j = 5J/J and d = 8D/J, we have 


= 8Q/Q 
where 


and 
R2 
C2 = sin? «— 4). 


jJ and d are small compared with J so that the fact that Cy and C2 vary by a 
few per cent during the life of the satellite will not be too serious and accordingly 


as a first approximation we may take 7 to be 7 042 km and « to be 65-28° in calculat- 
ing these coefficients*; hence 


Ci = 02511 
C2 


O°3119 
Merson & King-Hele’s result is that ()(theoretical)/Q(observed) equalled 
1°0068 + 0:0003 (standard deviation). The uncertainty of the equation for j and 
d arises largely from the uncertainty in « which is 0-02° in 65°. 
Hence 
10008 j —0-312d = (— 6-8 +0°3) x 1073, (1) 


Jeffreys’s adjustment of terrestrial data gives 


= (0+ x 1073, (2) 


A further equation for j and d can be obtained from Jeffreys’s (1942) study of 
the Earth’s gravitational field. The value of D given above is that used in the 
International Gravity Formula; it is not derived from observation but calculated 
on the hypothesis that the geoid is an exact spheroid (on this assumption the potential 
when expressed in spheroidal coordinates contains no fourth harmonic) and has 

* These are the values about half-way through the life of Sputnik 2 (Cornford 1958). 
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the value 4e(7e—5m), where e and m have their usual significance. The coefficient, 
y4, Of P, (cos 6) in the expansion of surface gravity is then 


4 
—gee(e— 5m) 
35 


or since e = J+4m, 


35 4 

If there is in addition a real fourth harmonic in the potential with coefficient 
5D, the contribution to y, is 8g,5D/7. 

Thus we find that 
8 8 
— = —(J-—2m)s] + -8D 
Se 35 7 
or 


Jeffreys’s harmonic analysis of free-air gravity gives 5y4/g, = (—0°7+ 6-4) x 
x 10-6 and hence 


I oy4 
= — = (—0°43 + x 1073. 
(—0°43 + 3°9) 

The uncertainty quoted for y, involves the neglect of the non-diagonal terms 
in the normal equations for the coefficients of the various harmonics and is there- 
fore a minimum estimate. Jeffreys found that the total variance of the nine fourth 
harmonic terms might be about 16mgal?; on this basis y, could be as much 
as 8-5mgal so that a standard deviation of 6-4mgal seems to cover the possible 
range of 4. 

With m = 3:450 x10-%, the observation equation is 


8 
1*23j x 1073 + = (—0°43 +3°9) x 1073 
or with quite sufficient accuracy, 
d= (—0°38 3°5) x 1073, (3) 


Two least squares solutions of equations (1), (2) and (3) have been made, one 
with weights inversely proportional to the squares of the uncertainties,* the other 
with equal weights. The results are given in Table 1. 

The unweighted solution gives a residual of 3-0 x 10~ or ten times the standard 
deviation for equation (1) and residuals comparable with and less than the standard 
deviation for (2) and (3) respectively. The weighted solution gives a very small 
residual for (1), one of 2-2 times the standard deviation for (2) and one o-g times 
the standard deviation for (3). In view of the few degrees of freedom on which 
all uncertainties are based, the weighted solution reconciles the data reasonably 


* The weights are 100 for equation (1), 2 for equation (2) and 1 for equation (3). 
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well. The uncertainties are based on 1 degree of freedom only and so are not very 
reliable. Nonetheless it is useful to compare the value of e with that derived from 
the constant of precession and from gravity data. The former method gives 
e~1 = 297°34 with a standard deviation of 0-28 or 0-55 depending on whether the 
Earth is assumed to be in a hydrostatic state or not (Jeffreys, 1948). The difference 
between the precessional value for e~! and the satellite value is 0-go, equal to 
twice its standard deviation on the hydrostatic hypothesis. There is thus some 
evidence of a real departure from the hydrostatic state. 


Table 1 
Solutions of observation equations 
Weighted Unweighted 

j (—5-74 + 2-6) x 107% —3:01 X 107% 
d + 38-2) x107% +2:63 x 1073 
—9°4 —4'92 X 1078 
8D x1078 —4°3 x 1078 

(1-628 + 0-004) x 1073 1-632 x 107% 
ew! 298-24 + 0°36 297°88 


The uncertainties are standard deviations. 


Gravity data alone give values of e™1 which vary according to the treatment of 
the possible longitude terms but are up to 2:0 less than the satellite value. It 
is important to emphasize that Jeffreys’s analysis of the Earth’s gravity field, 
from which the coefficient y, is obtained, is out of date. It depends on data 
available in 1940 and since then a great deal more material has accumulated, 
especially in the Pacific Ocean. In view of the precision attainable with satellite 
observations, it is very desirable to analyse this material using Jeffreys’s statistical 
procedure and so obtain estimates of J and D with uncertainties more comparable 
with those of satellite observations. The importance of this is the greater in that 
observations of the Moon and the calculation of J from the precessional constant 
and the Earth’s internal density distribution cannot be expected to give much 
greater precision than at present. 

This note has been prepared for the Satellite Study Group of the Geodesy 
Sub-Committee of the National Committee for Geodesy and Geophysics, and is 
published with the permission of the Director of the National Physical Laboratory. 


Standards Division, 
National Physical Laboratory, 
Teddington. 
1958 September 22. 
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Letter to the Editors 


The Eccentric Geomagnetic Dipole 
W. D. Parkinson and J. Cleary 


The Gauss coefficients of a spherical harmonic analysis of the Earth’s main 
magnetic field have been published recently (Finch and Leaton 1957) for the epoch 
1955.0. The authors state the positions of the poles of the centred dipole (usually 
called the geomagnetic poles) but not those of the eccentric dipole. The positions 
of the latter were required here recently, and are presented in case they may be of 
interest to others. Phenomena not involving the surface of the Earth, such as the 
deflection of auroral or cosmic ray particles outside the atmosphere, are controlled 
by the eccentric rather than the centred dipole. 

Formulae for deriving the position of the eccentric dipole from the Gauss 
coefficients have been given by Bartels (1936). On substituting the 1955.0 values 
into these formulae, we obtain the following results: the eccentric dipole is dis- 
placed by 0-0685 Earth radii (about 436km) from the centre towards a point at 
latitude 15°-6 N, longitude 150°-9 E (just east of the Marianas). This is 6°-6, 
or 730km north of the geomagnetic equator. The poles of the eccentric dipole 
(i.e. the points where its axis cuts the surface of the Earth) are at 81°-o N, 84°-7 W 
and at 75°-0 S, 120°-4 E (i.e. in Ellesmere Is. and in Wilkes Land). The axis of 
the eccentric dipole is not vertical at these points, but is inclined at 3°-9 to the 
vertical in the direction of the corresponding geomagnetic pole, i.e. 126°-3 E of N 
for the north pole and 150°-9 west of north for the south pole. 

The field of the eccentric dipole is, of course, parallel to its axis, and so it also 
is inclined to the vertical at the poles of the eccentric dipole. The eccentric dipole 
field is vertical at two points which are further from the geomagnetic poles. The 
positions of these points are 82°-4 N, 137°:3 W and 67°-9 S, 130°°6 E. The latter 
is close to the observed dip pole near Terre Adelie, but the former is about 1 200km 
from the observed northern dip pole which is near Prince of Wales Is. (Hope 1957). 


Bureau of Mineral Resources, 
Geology and Geophysics, 
Department of National Development, 
Melbourne, Australia. 


1958 July 29. 
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Notes on the Progress of Geophysics 


Geophysical Studies in Polar Regions 
The Antarctic Ice Sheet 
G. de Q. Robin 


Introduction 


The concept of the International Geophysical Year has done more than 
capture the imagination of scientists. The problems of placing men and materials 
at the most inaccessible spots on Earth have produced an enthusiastic response 
from transport engineers and explorers while adequate financial support has 
enabled them to achieve their objectives. Still more remarkable are the recent 
voyages of nuclear powered submarines under the Arctic pack ice, which open up 
new and more rapid techniques for the exploration of this oceanic basin. It is 
therefore not surprising that 1957-58 has seen the available geophysical data on 
the polar regions increase by an order of magnitude. 

Although the major disciplines of the IGY are those most directly linked with 
solar radiation in its various forms, certain studies of the structure of the Earth 
were added as the programme was planned, and this review will be limited to some 
of the latter investigations. Most of this work is still in progress, but brief sum- 
maries have been published which are sufficient to outline the principal results 
to date (IGY Bulletins 1957-8, Hatherton 1958, IGY Glaciological Rep. 
1958, Shoumsky 1958). 

The main advances in the study of the geophysics of the Earth’s crust in polar 
regions have come from the application of existing techniques of geophysical 
exploration. Such techniques have proved particularly suitable for studying the 
ice masses of the polar regions. Field studies have followed the laboratory and 
theoretical work of the past decade on the flow of ice, with the result that glaciology 
is now recognized as a subject for geophysical study which, in addition to its own 
intrinsic interest, may help provide solutions to broad geological problems on the 
flow and deformation of rocks. 

The basic data on the depth and other properties of the polar ice sheets have 
been gathered by well equipped scientific traverse parties. Seismic reflection shoot- 
ing methods have provided the main key to ice depth determinations while gravi- 
tational and magnetic measurements are used to fill in details between the main 
seismic stations. 


Seismic shooting methods 


The contrast in density and elastic wave velocity between ice and rock is so 
favourable for seismic reflection shooting that rapid exploratory surveys making 
widely separated spot measurements have proved successful in most cases. Occa- 
sional areas are found where the reflecting properties of the glacier bed are poor, 
as in North Eastern Greenland, and here one has to rely on gravity measure- 
ments to a much greater extent (Hamilton 1956). However, frequent double or 
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triple reflections from the glacier bed confirm that the ice-rock interface usually 
forms the good reflector which one would expect on theoretical grounds and also 
show that attenuation of elastic waves in ice is relatively small at the frequencies of 
the order of 1ooc/s normally recorded. Ice thickness measurements made by 
reflection shooting now vary between 100m and 4 200m and the accuracy of depth 
measurement is normally within two per cent. The technique is therefore adequate 
but any simplification of the complicated equipment and lengthy procedures in- 
volved would be welcome. Unfortunately the complications increase as one climbs 
onto the high inland plateau of Antarctica. Whereas one kilogramme of explosive 
on or near the surface is sufficient to produce adequate reflections at lower altitudes, 
on the cold inland plateau such explosions produce large and prolonged surface 
waves which tend to block out weak echoes. (Robin 1958). The exact nature of 
these surface waves is not yet clear, although they are obviously connected with a 
very high velocity gradient close to the surface of the cold inland ice. A wave- 
guide effect appears to account for some of their characteristics, but does not explain 
the sudden increase of surface wave energy found when mean ice temperature 
falls below approximately —30°C. The only reliable way to avoid this difficulty 
is to increase the shot depth to 20 or 30m so that very little energy is captured by 
the surface “‘waveguide”. This has been done with conspicuous success by the 
Australian traverse team in the Antarctic, while other traverse teams on the 
Antarctic plateau with less effective drilling equipment have encountered varying 
degrees of difficulty. 

On the floating ice shelves around the Antarctic continent, the theoretical 
expectation of a reflection coefficient of energy of P waves of around 0-16 has not 
generally been fulfilled in practice, although this figure has been found to apply 
on some occasions (Robin 1958). The most reliable technique for this problem 
has been to employ a combined refraction-reflection technique using the reflection 
from the sea bed. 

Refraction shooting techniques using spreads up to 20km have been used 
successfully by the Russians for measuring ice thickness and wave velocities in 
the underlying rock (Shoumsky 1958). However in spite of the extra data obtained 
this technique is slow and uses too much explosive for most exploratory surveys 
of Antarctic ice. On the small scale, refraction shooting over spreads up to a few 
hundred metres have been used extensively to determine the velocity gradient, 
and hence the related density variations in the top 50 to 100m of the ice sheet 
where snow is changed gradually into glacier ice. (IGY Glaciological Report 
1958). 

The absence of any discontinuities in velocity in this upper zone and the purity 
of the ice—air mixture make such polar ice sheets a useful testing ground, inter- 
mediate in size between laboratory and earthquake studies, on which experiments 
on elastic wave problems under accurately known conditions can be carried out. 
It is to be hoped that the use of three-component seismographs and absolute meas- 
urements of energies under these conditions will provide information which will 
help to consolidate the bridge between certain theoretical and experimental seis- 
mological studies, particularly in the case of amplitude problems. 


Seismology 


In the field of earthquake seismology, a scattered network of stations on the 
Antarctic continent has been established. Stations have been operated success- 
fully on the thick inland ice (e.g. South Pole) and on a floating ice shelf (Halley 
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Bay). No major earthquake activity within the continent has been reported 
(Hatherton 1958), although smaller shocks have been recorded at Halley Bay 
(Weddel Sea) and at Scott Base (Ross Sea), the latter including activity associated 
with the active volcano Mt. Erebus. Many shocks from elsewhere in the southern 
hemisphere were recorded, particularly from the ridge in the Indian Ocean. 
An interesting feature of microseismic activity was the very low level to which this 
fell during the winter months at Scott Base, and to some extent at D’Urville Base 
further north. The explanation may be found to be an effect of the pack ice 
which covers greatly increased areas of the surrounding ocean during winter 
months. Another promising field of investigation from which results are awaited 
are studies of surface waves since these may give information on the mean thickness 
of the Antarctic ice sheet and thus form a check on the representative character 
of the many spot determinations of ice thickness. 


Gravitational studies 

Apart from a very few basic measurements of gravity which have been checked 
at points on the Antarctic coast line, the main activity has been to use gravity meters 
for interpolation between seismic soundings of ice thickness. In spite of the diffi- 
culty of measuring elevations by aneroid barometers, the method has proved most 
successful for providing a rapid method of ice depth measurement, even though 
its accuracy is less than that of seismic soundings. As a change of 1 mgal in the 
value of gravity is caused by a change of ice thickness of around 14 m, an elevation 
error of 20m will cause an error in ice thickness of about 100m. Local gravita- 
tional anomalies arising from changing densities of rocks may well produce larger 
errors. Elevation measurements should be sufficient for interpolation of ice thick- 
nesses between seismic soundings, but it is difficult to eliminate the possibility of 
systematic errors occurring over the vast distances between the periphery and 
central regions of the Antarctic continent when measuring elevations by baro- 
metric pressure. For this reason it is difficult at this stage to interpret the reported 
negative free air anomalies of the order of one hundred milligals on the Antarctic 
plateau. Even one major levelling traverse using optical methods, of the type 
carried out in Greenland by the British North Greenland Expedition (Paterson 
& Slessor 1956), would be of considerable value in the Antarctic. In north Green- 
land free air anomalies were generally less than 30mgal over the central portion 
of the ice sheet (Hamilton & others 1956) thus showing that isostatic compensation 
for the weight of the ice sheet had taken place. As the establishment of such 
isostatic equilibrium is believed to take thousands of years, this indicates that this 
portion of the Greenland ice sheet has been approximately stable for many cen- 
turies. Apart from its considerable scientific significance, similar data for Antarctica 
would provide a tangible return for those administrators who were persuaded to 
provide financial support for Antarctic research because of the widespread damage 
which would be caused by sea level changes resulting from a decrease of Antarctic 
glaciation. 


Magnetic observations 


Although several traverse parties have taken regular measurements of magnetic 
intensity, particularly of the vertical component, no preliminary reports on their 
value have been seen. It appears likely that apart from general confirmation of 
gravitational anomalies, they will be both less accurate for ice thickness inter- 
polation and of less general significance than gravitational observations. 
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350 G. de Q. Robin 
General glaciological problems 


In order to understand the mass and thermal balance of the Antarctic ice sheet 
a vast amount of data is being collected by many traverse parties. The net annual 
surface accumulation of snow has been found to vary from 40cm of water equiva- 
lent in some coastal areas to 8 cm or less on the high inland plateau. Ice streams 
draining the inland area have been found to move up to 3m per day, (IGY Glacio- 
logical Rep. 1958, Shoumsky 1958), but owing to the difficulties of access satis- 
factory estimates of the rate of discharge of the largest ice streams have still to be 
made. Until this is done this method of studying the mass balance will remain 
ineffective. However other evidence now available, such as the domed shapes 
of the ice sheet, the scarcity of freshly exposed rock surfaces and the small values 
of the world wide changes in mean sea level, all point to the existence of an approxi- 
mate mass balance of the Antarctic ice sheet. 

Some headway is also being made in understanding the temperature distribu- 
tion in the Antarctic ice sheet thanks to deep drilling projects carried out by the 
U.S.A. and the U.S.S.R. Both boreholes showed decreasing temperatures with 
depth near the surface, but the Russian results near Mirny showed that this 
effect was reversed at around 200m depth so that the basal ice temperature was 
close to the pressure melting point at a depth of just over soom. They have used 
the method of hydro-analogies to study the temperature distribution in the ice, 
and have obtained reasonable agreement with observations (Shoumsky 1958). 
This should provide an improved method of analysis, although the advection of 
cold ice remains the primary explanation of these negative gradients (Robin 
1955): 

A continuous trend of climatic change may also affect observed temperature 
gradients considerably, but in order to eliminate the effect of horizontal advection 
of cold ice, boreholes to study this problem must be sited at a centre of outflow of 
ice. So far no boreholes have been shown to satisfy this criterion, hence this form 
of a direct physical record of climatic change has still to be obtained. 

The U.S.A. boring at Byrd base was also notable for a 98 per cent core recovery. 
Half the core has been returned to the U.S.A. for detailed laboratory studies which 
should prove very valuable, as it will provide a well dated sediment of the past 


700 years for studying problems such as the rate of fall of micrometeorites and other 
particles. 


Physiography 

At the beginning of 1956 an effective 400 miles or so of ice thickness profile 
had been obtained from Antarctica. At the time of writing over six thousand miles 
of profiles have been obtained and further results are anticipated. It now seems 
likely that the mean thickness of the Antarctic ice sheet is around 2000m. While 
from a geological point of view it has always been clear that the Antarctic block is 
of continental proportions, some results have again raised the physiographical 
question of how much of the rock surface lies above sea level. It is now clear that 
that boundary of East Antarctica which faces West Antarctica is above sea level 
for at least 2 000 miles, a figure of continental proportions, but two major problems 
remain. The first is to prove or disprove the existence of an ice filled trench 
between East and West Antarctica. The evidence to date tends to confirm the 
existence of this discontinuity which has been postulated by geologists for many 
years. The second major problem is to find out whether or not the huge ice mass 
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on East Antarctica has depressed the rock floor below sea level to form a huge 
basin-shaped depression. It seems reasonable to hope that the U.S.A. and U.S.S.R. 
will provide the answers to these problems in the 1958-59 field season. 


Crustal geophysics 

So far no seismic shooting or seismological results have been reported which 
indicate the depth of the Mohorovicié discontinuity beneath the Antarctic ice 
sheet. The study of this problem may provide crucial evidence on whether or not 
this discontinuity is due to a phase change. If there is a change of material, a 
superimposed ice mass 3 km thick would only depress the discontinuity by around 
1km, which is not significant. On the other hand much larger changes may be 
possible if the discontinuity is a function of temperature and pressure as would be 
the case if the discontinuity is the product of a phase change. Both theoretical and 
experimental methods are needed for the solution ot this problem. 


Conclusion 


In concluding this review, one must pay tribute to the vast number of scientists 
and technicians who have put up with the extreme discomforts and hazards of the 
polar regions to gather this vast collection of data which is advancing our knowledge 
of the polar regions at an unprecedented rate. 


Scott Polar Research Institute 
Cambridge. 
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The 83rd Colloquium of the French National Centre of 


Scientific Research (Topography and Ecology of Oceanic 
Deeps) held at Nice, 1958 May 5-14 


Representatives from Norway, Sweden, Denmark, Holland, U.S.A., Italy, 
Spain, North Africa and Great Britain attended this meeting, and a most useful 
exchange of ideas was possible, both between oceanographers from different 
countries and between specialists in subjects whose domains overlap. Not 
unnaturally, the Western Mediterranean basin was the subject of much discussion, 
and the combination of geological and geophysical work has produced some in- 
teresting ideas on the origin of the Mediterranean. It was an opportune moment 
for clarifying various viewpoints on Mediterranean structure, since not only is 
there considerable drilling activity in the area, but also the U.S. research ships 
Vema and Atlantis have been making a series of seismic measurements in the Medi- 
terranean. 

Full transcripts in French of all the papers have been made by the Colloques 
Internationaux du Centre National de la Recherche Scientifique (13, Quai Anatole- 
France, Paris 7°). 

The colloquium at Nice will probably be remembered as marking the begin- 
ning of the bathyscaphe era in submarine geology. Slow development has been 
going on since the end of the war by the French Navy, the Belgians and the school 
led by Professor Picard. Funds for it have always been short, but successful 
dives have been made in the deep part of the Mediterranean and the difficulties 
and hazards of this type of work have been overcome. The U.S. Navy is now in- 
terested in bathyscaphes and is acquiring some for their own experimental work. 
There is little doubt that a rapid growth in bathyscaphe exploration of the sea-bed 
will soon follow. 

The bathyscaphe enables oceanographers to have a direct look at what is beneath 
the sea surface and, in particular, what is on the sea bed. It should be possible 
to study deep-sea marine life with the help of the bathyscaphe. Knowledge of the 
animals that inhabit this domain is at present restricted to those bodies that can be 
brought to the surface in dredges and nets, and the habitat and behaviour of the 
animals is very much a matter of speculation in the dark. Photographs of the sea- 
bed have enabled marine geologists to formulate some ideas of the form of the 
deep-sea floor, and in some places ripple marks show the existence of currents in 
the water. However, a series of single photographic exposures, with a limited 
field of view, are not so useful as a continuous look from a vessel moving slowly 
at a short distance from the sea-bed. The shape of the sea-bed is known roughly 
from echo-sounders, but the detail is blurred because the echo-sounder is normally 
three miles above the sea-bed and records the depth to the nearest point, which 
may not necessarily be vertically beneath the instrument. It will be possible for 
the bathyscaphes to recognise whether topography that looks on the echo-sounder 
like hills and valleys is in fact similar to such features on land. 

The bathyscaphe can obviously be of use in sample collecting. One is never 
sure that the corer goes into a typical part of the sea-floor, and much time and 
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effort could be saved by visual selection of the site. Current measurements could 
be made from the bathyscaphe at the sea-bed, and any movement of the sea-bed 
caused by storms on the sea-surface will be observed. 

Professor L. Glangeaud described the Mediterranean as being neither an epi- 
continental sea like the Japan or China sea areas, nor an ocean like the Atlantic 
and Pacific with typical 6-7km/s material close to the surface of the sea floor. 
From a geological study of the margins of the Mediterranean Professor Glangeaud 
has produced a mechanism of evolution for the Mediterranean which explains 
many of the geophysical and topographical data. The Mediterranean originally 
could have been a true deep ocean separating two continental sialic blocks— 
Europe and Africa. Normal denudation and sedimentation increased the size of 
these continents and, at the end of the Carboniferous period, the growth of the 
European and African Hercynian chains completely filled the gap separating 
Europe and Africa. Thus during the upper Permian and the lower Trias, an 
emerged continent stretched from England to the Sahara. There remained 
perhaps a few elements of the primary Paleotethys in the eastern Mediterranean, 
but there is no evidence to prove it. The sea reappeared from middle Triassic 
times and there is evidence which suggests that it did so because of a drifting 
apart of the European and African blocks. While such drifting is only postulated 
to meet the demands of surface geology, it is hoped that palaeomagnetic measure- 
ments may one day provide positive evidence. There is an increasing body of 
evidence that favours drifting of land masses; the mechanisms for such drift will 
no doubt be forthcoming when it has been demonstrated that they have, in fact, 
taken place. The latest evidence shows that rift valleys in the North East of 
Africa are formed by tearing apart the Earth’s crust rather than by compressional 
forces and consequent pushing down of a block of material. If such drifting 
apart was a feature in Mediterranean evolution, it was followed by a phase in which 
the new sea again filled up and was raised in the Oligocene for a short period. 
From the time of the Miocene, sinking occurred and the sea covered the debris 
of the Alpine chain from the Carpathians in the north to the Saharian in the south 
—in much the same way as the Triassic sea had covered the debris of the Hercynian 
chain. 

If this picture is substantially correct, the seismic measurements should show 
the Moho to be at a depth of ten miles or less below the sea-surface. The few re- 
fraction stations reported by Gaskell and Swallow in 1953 showed the presence of 
at least 2}km of rock in which the compressional wave velocity was 44 km/s. 
The energy propagation was bad and rough weather caused bad background noise 
during the observations and no 6-7km/s refracting layer was detected. Uncon- 
firmed reports from the U.S. Mediterranean expedition of 1958 suggest that the 
same bad propagation was experienced and that the structure of the deep Medi- 
terranean basins is still undetermined. No doubt there will be some depression 
of the original ocean floor, due to the abnormal amount of material deposited 
since the Mediterranean was a true ocean. In fact, although Professor Glangeaud 
asserts the contrary is true, it is possible that the Mediterranean is similar to the 
Japan sea region, where several kilometres of rock in which the compressional 
wave velocity is about 5-8km/s are found below the sea-bed instead of the usual 
6-7km/s found in other oceans. 

Professor Kuenen, on the other hand, was of the opinion that the Mediterrane in 
is an example of turning continental into oceanic type by some form of melting or 
dissolving of the granitic layer, together with the dropping of a whole block. 
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There certainly needs to have been a large mass of material to provide the sediments 
of the Alps and the Sahara. 

The topography of the deep oceans was discussed by Dr Menard of the Scripps 
Institution of Oceanography. Menard has calculated that there are about 5—10 000 
atolls, volcanos, guyots and sea-mounts in the Pacific. This figure is estimated 
from the value of 800 already discovered and the fact that only 10 per cent of the 
Pacific has been explored; from the length of track of ships with echo-sounders 
and the value for number of features per mile of track; and thirdly by taking a 
well-explored area and extrapolating. This large number of sea-bed features 
means that the Pacific consists of a series of small hills protruding from a flat 
base. Around the larger features are smooth areas that Menard has called archi- 
pelagic aprons—these are large compared to their parent volcanos and therefore 
are probably formed by erupted volcanic material rather than by sediment derived 
from the erosion of the volcanos. On this picture, the area of sea bed smoothed by 
volcanic outpourings is as great as the flat areas produced by submergence of 
sea-bed features by sediments moved by turbidity currents. The same topo- 
graphical features probably occur in other oceans, but so far only a small line of 
flat topped sea-mounts has been found in the Atlantic. The deep-ocean seismic 
results summarized by Gaskell suggest that there is a difference between the 
Atlantic and the Pacific. The Atlantic—at least over the greater part of its area— 
can support volcanic features that grow from the sea-bed and it is possible that this 
is due to a hardening of the rock that comprises the sea-bed. It is significant that 
the places where hardening of the rock is needed to explain the non-sinking of 
volcanic islands are just those where continental drift has been postulated. It could 
be that movement of continental areas leaves an altered 6-7 km/s layer (or sima 
layer) which is more rigid than the original ocean floor. On this hypothesis, of 
course, the continental drift is irreversible because it can only take place into the 
slightly plastic primaeval ocean floors. 

Dr Heezen of Woods’ Hole showed his new charts of the underwater topography 
of the Atlantic Ocean. These charts are made so that the different sea-bed topo- 
graphical features—mountain ranges, isolated volcanos, flat abyssal plains, etc.— 
stand out clearly. They may be obtained from The Geological Society of America, 
419 West 117 Street, New York 27, N.Y. Part I of Special Paper 65, Physiographic 
Diagram of the North Atlantic, Sheet 1, by Bruce C. Heezen and Marie Tharp is 
issued, price $1.50 per copy. Part 2, the text, by Heezen, Tharp and Ewing, 
consists of 150 manuscript pages, 49 figures and 30 plates, and is in press. Dr 
Heezen suggested that the evidence of such features as the rift in the Mid-Atlantic 
Ridge, and the recent evidence that the African rift valleys are due to tension rather 
than to compression, indicate that the Earth is expanding, and therefore a warming- 
up, not a cooling-down, is taking place. Warming-up is postulated by those who 
would build the Earth by collection of material from outside the Earth. A recent 
paper in the Geophysical Journal by Dr Lubimova (1958) maintains that the com- 
bination of radioactivity and poor conductivity of the rocks that constitute the 
Earth must mean that the inside is getting hotter—although the outside crust could 
have cooled in the past. This could well explain Heezen’s rift zones. 

Sediments in the deep sea were the subject of Dr Wiseman’s paper. Dr Wiseman 
stressed the need for great care in taking of deep sea cores. Apart from the technical 
point of disturbing the layering of the core by unsuitable design of coring device, 
the actual location of the core is of paramount importance. For example, cores 
from abyssal plains and from continental rises are unsuitable for studying past 
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temperature changes in the sea because the sediments will have been subjected to 
slumping and to erosion. Movement of sediment by such burrowing animals as 
holothurians or by ‘sand volcanos’ must also be avoided. Dr Wiseman suggests 
therefore that cores should be taken from the top of an isolated knoll. 

Professor K. O. Emery of the University of Southern California described his 
work on the sediments of the continental borderland of S. California. The rate 
of sedimentation is measured by making C analyses on core-samples and ranges 
from 20 mg/cm?yr at the bottom of the basin to 60 mg/cm?yr nearer land. A rate 
of 50mg/cm?yr corresponds to about 1 mm thickness of deposit per year at the 
sea-bed, so that these figures are about ten times the values found in the deep 
oceans. Some care must be taken in the carbon age determination because the 
age at the top of the core is found to be 2000 years, and this value must be sub- 
tracted from the age of sections of the core lower down in order to obtain the cor- 
rect deposition rate. 

Dr Landergren proposed the use of the measurement of the boron content of 
sediments in order to determine salinity conditions when sediments were deposited. 

A visit was paid to the National Museum of Oceanography at Monaco. This 
museum has an aquarium and a fine collection of old oceanographic devices. In 
addition to modernizing the lay-out of the exhibition, the Director, Commandant 
Cousteau, proposes to follow the wishes of Prince Albert of Monaco, the founder 
of the museum, by doing as much deep-sea oceanography as possible. 

The museum hopes to develop techniques which will enable difficult experi- 
ments to be carried out at sea. The underwater diving with the aqualung has proved 
of great value in studying underwater geology, and Commandant Cousteau described 
the work he is now doing on the use of nylon cables for dragging dredges and for 
anchoring ships in deep water. 

Two of the three directors of the International Hydrographic Bureau presented 
papers during an afternoon visit. Director Viglieri described the progress of the 
latest revision of the Monaco Hydrographic maps and spoke of the co-operation 
that has been arranged with the various national hydrographic services. Director 
Damiani discussed the present state of agreement and disagreement on the subject 
of nomenclature of ocean bottom features. 

The Colloquium closed with the official inauguration of the Oceanographic 
Station at Villefranche-sur-Mer (Director: Professor J. Bourcart). 
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Book Reviews 


Principles of Geodynamics 
Adrian E. Scheidegger 
xi + 280 pp. (Springer-Verlag) 


Dr Scheidegger remarks in his preface “It seems, therefore, that the time is 
ripe for an evaluation of existing ideas (on geodynamics) in the light of presently 
available facts’. This is undoubtedly so, and it is time for a critical review of 
facts and theories about such topics as the origin of continents and oceans, the 
formation of geosynclines and mountains, and earthquake and volcanic activity. 
Unfortunately this book just misses being the review that is needed. The title 
itself is rather misleading. The successive chapters are summaries of physio- 
graphical, geological and geophysical facts concerning the Earth; of the mechanics 
of deformation and theories concerning the effects of the rotation of the Earth and 
of ideas that have been put forward about the formation and growth of continents 
and oceans, orogenesis, faulting and folding and miscellaneous minor features. 
It is, however, disappointing that there is little discussion of principles, in the 
sense of conditions to which all theories should conform—a sort of geophysical 
“fence” to bound future arguments. It might be thought that these conditions are 
obvious, but in view of the evident difficulties of the subject, stemming from the 
fact that we can observe only the state of the surface of the Earth at the present 
moment of time, it is surely worthwhile to attempt to make these principles 
explicit. Almost everyone, for instance, tacitly assumes that the behaviour of the 
Earth is governed by the laws of classical physics, although Egyed has argued 
that some phenomena are the consequence of changes in the ‘‘constants” of physics 
that, it has been suggested, are entailed by general relativity. (Professor Wilkinson 
has adduced evidence from pleochroic haloes to argue that changes in fundamental 
atomic constants are extremely small (Wilkinson 1958).) 

Another principle to which most theories seem to conform is that the state of 
the Earth in the past should be assumed to be much as it is today. Geologists 
especially have been very successful in interpreting their data on this basis, a 
success which may have obscured the fact that classical geology has dealt only with 
the last third or quarter of the Earth’s history and that the origin of the continents, 
for instance, seems to be early in the Earth’s history. In one respect at least it is 
untrue that the present state is representative of the past, for if the usual laws of 
physics are supposed to apply, then the present stock of radioactive matter is less 
than the original. Estimates of the importance of this factor depend on the 
assumed initial temperature distribution in the Earth, since this determines the 
proportion of the present efflux of heat that is due to radioactivity. Our views on 
the heat energy available for geodynamic processes at any time in the Earth’s 
history must therefore depend on our views about the origin of the Earth. In this 
connection, one would welcome a discussion of the sources of available energy. 
It usually seems to be assumed that these are either stored heat from the Earth’s 
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initial state or energy stored in radioactive nuclides. There is also the kinetic 
energy of the Earth-Moon system and, if there are phase changes within the Earth, 
energy would have been released from low-pressure forms of minerals as the Earth 
condensed (whether from gas or from solid particles). It would, incidentally, be 
useful to know the time constant for such phase changes. 

The only section of the book in which principles in this sense are discussed is 
the chapter on the mechanics of deformation. The importance of a very careful 
study of this subject cannot be overemphasized, and sophisticated summaries are 
given of the theories of finite strain, of infinitesimal elastic deformation and of 
some ideas of plasticity and hydrodynamics. Unfortunately the parts of the subject 
with the most sophisticated development are those that are of least use in geodyna- 
mics and the balance of this chapter, although it reflects the present state of know- 
ledge of the deformation of solids, does not reflect the relative importance of the 
parts of that subject to geodynamics. The primary difficulty here, as in most 
aspects of geodynamics, is the gross discrepancy between the time scale of natural 
events and that of anything we can observe. Accordingly, although various effects 
may be identified as resulting from certain types of rheological behaviour, an enor- 
mous extrapolation is involved and it is arguable that it may in future be possible 
to establish the rheological behaviour theoretically from the theory of solids with 
more certainty than it can be inferred from experimental observations—this may 
even now be the case in the analogous problem of estimating the equation of state 
of matter in the core of the Earth. Dr Scheidegger gives us a useful summary of 
the evidence for the types of rheological behaviour in the Earth’s crust for different 
time scales but part of this is now superseded since Sir Harold Jeffreys has shown 
that the law propounded by Lomnitz fits various observations rather well. But not 
only is the physical behaviour of matter over long intervals very poorly known; 
in addition the mathematical theory of the motion of bodies following rheological 
laws that may hold for the Earth has important gaps. Dr Scheidegger mentions 
the most serious of these gaps—the study of rheological instability—and remarks 
that mountain ranges are probably the manifestations of lines of instability. 
Bijlaard has considered instability of a plastic material in tension but this is 
clearly irrelevant to the problem of orogenesis, where it seems that instability in 
compression occurs. It is likely that allowance must be made for the weakening of 
the crust due to heating by the blanketing of sediments once a geosyncline has 
started to form. Another example of instability is probably the flow of rocks pro- 
ducing the folds characteristic of mountains. The patterns of these folds are very 
reminiscent of instability in hydrodynamic flow or magneto-hydrodynamic pheno- 
mena where a displacement gives rise to forces tending to increase the curvature 
of the flow line or line of force. 

In view of Dr Scheidegger’s aim of confronting the various theories of geo- 
dynamical problems with the available facts, one would have wished for a more 
extended critical discussion of these facts instead of the rather condensed summary 
that is given. This, to my mind, is the main weakness of the book. It is very useful 
to have the theories most competently summarized, but one feels that with a fuller 
discussion of the facts and their significance this book would have begun a new 
advance in these studies. 

One point that might have been discussed in more detail with quantitative 
data is the significance of isostatic balance and the departures from it, particularly 
so far as inferences can be drawn about the strength at different levels in the Earth. 
‘The very similar rates of flow of heat from continents and oceans is another 
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important fact, while the different depths of the Mohorovici¢ discontinuity below 
continents and oceans and the different chemical constitutions of the crust deserve 
more than mere statement. But attention should also be called to the differences 
between one ocean and another, such as the clear distinction between the topo- 
graphy of the floor of the Atlantic and Pacific oceans. The mass of detail we have 
about the continents may obscure equally great distinctions between continents, a 
striking instance being the great mountainous areas of Central Asia which have 
no parallel elsewhere and which find no explanation in any of the usual theories of 
mountain building. It is not Dr Scheidegger’s fault that his section on palaeo- 
magnetic data will certainly soon be out of date, for this is a rapidly developing 
subject; the data so far come entirely from the continents but in view of the 
inference of continental movement that has apparently to be drawn it would be 
worth considerable effort to attempt to obtain similar data from oceanic regions. 

In considering the origin of mountains, Dr Scheidegger introduces a relation 
between the cross-sectional area of a mountain range and the apparent shortening 
suffered by the rocks. The relevance of this apparent shortening, inferred from the 
folding of the visible rocks, requires more critical discussion. One may certainly 
distinguish between the deformation of the crust and that of the much weaker 
geosynclinal deposits above it and it has already been suggested that the pattern 
of the latter is comparable with patterns of unstable flow. Geologists have, in fact, 
sometimes suggested that nappes result from the sediments sliding down the slopes 
of a rising crust. If these ideas have substance, then the apparent shortening is 
irrelevant to the main problem, that of the deformation of the crust. ‘This might 
profitably be considered in terms of energy balance. The excess potential energy of 
a mountain may be estimated from its topography and gravity anomalies and the 
shortening may be estimated from the supposition that the cross-sectional area is 
unchanged in the deformation. Then the maximum work done is the product of 
the shortening by the strength of the undeformed crust and this must equal the 
potential energy of the mountain range plus any dissipation in plastic deformation. 
The latter is the most difficult to estimate and more theoretical work on losses in 
plastic deformation would be required; this, by the way, might be a more profitable 
field for model studies than investigation of the geometry of deformation. It can 
however be estimated that the latent heat of a large batholith is only a small fraction 
of the potential energy of a typical mountain range so that the formation of batho- 
liths may, so far as energy considerations go, be regarded as incidental to the 
formation of a mountain. This is a very simple-minded study proposal and a 
considerable difficulty in putting it into effect would be deciding to what depth 
the deformation can be considered to extend. 

One may think of the crustal boundary as a chemical or phase boundary, a 
density boundary or seismic boundary, a stress boundary or a strain boundary. 
The seismic boundary is the only one that is well established and it is more than 
a little naive to suppose that all the others coincide with it and form a well-greased 
surface over which the crust may slide. The association of deep focus earthquakes 
with geosynclines shows that there is no firm stress or strain boundary at the 
Mohorovicic discontinuity and a deep seismic boundary is not always found 
where indicated by gravity anomalies. 

Of the many concepts of orogenesis, that of the “‘tectogene” as a buckled 
region in the crust should certainly be forgotten. It should never have been given 
currency as the idea of a region in which the inclination of the crust was so great 
that the engineering theory of elastic buckling applied. This probably arose from 
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misleading diagrams with exaggerated vertical scales, but in fact there is no evidence 
that the crust deforms in this way and Dr Worzel has now shown that gravity 
anomalies over island arcs indicate a thinning rather than a thickening of the 
crust. In considering the development of a geosyncline, it is probably wrong to 
think of sediments accumulating uniformly along the length of the geosyncline. No 
other geosyncline has been studied in such detail as the Welsh geosyncline and 
there Jones’s (1956) results show how at one period the sediments accumulated 
with their greatest thickness first in one part, then in another of the geosyncline, 
until the whole geosyncline was filled up. The thick sediments were almost 
entirely deposited in shallow water so that one has to think of the major depression 
of the crust occurring in different places at different times, with strong indications 
that the local rate of change of depression is too great to be accounted for by elastic 
flexure of the crust. 

Dr Scheidegger makes one point in his discussion of theories of orogenesis 
that is questionable. He says (p. 161) that metamorphism is associated with 
increase of density. If, however, the formation of granite is regarded as the final 
stage of metamorphism, as Dr Scheidegger seems to accept, then this is not true, 
for it is well known that granites are almost always less dense than the country 
rocks, even where the latter are scarcely more than highly compact sediments. 

Another question to which more attention might be given is classification of the 
types of deformation that are found. Mountains are fairly obviously the result of 
compression while it seems very likely that rift valleys are due to extension, and 
now current work on earthquakes round the Pacific suggests that the predominant 
motion is a rotation of the Pacific relative to the surrounding land. 

This is certainly a useful book bringing together all the current and past theories 
and confronting them with observation at a time when new discoveries, such as 
those that may be made when a borehole passes through the Mohorovici¢ discon- 
tinuity, may turn the subject upside down. A very valuable feature is the exten- 
sive reference to recent Canadian and Russian work and the book owes much to 
the author’s association with Toronto and Ottawa. While it is remarkably up-to- 
date, containing for example a number of references to papers read at the I.U.G.G. 
meetings at Toronto in September, 1957, the book also shows some surprising 
omissions from the literature, for instance the work of Jeffreys and Bland (1957) 
on convection in the mantle of the Earth. There are also indications that the author 
was not clear for whom he was writing, for some things, such as the principles and 
methods of seismology, are explained in a very elementary way while others, like 
the theory of finite deformation, receive an advanced treatment. 

The printing and diagrams are excellent and there are full author and subject 
indices, but the English reader may often be brought up sharply by the use of 
“lap” meaning “‘laplacian’’, while there are several at.noying minor misprints. 

There is no doubt that this book deserves careful study by everyone with a 
serious interest in its subject. 


A. H. Cook. 


References 


Wilkinson, D. H., 1958. Phil. Mag. (8) 3, 581. 

Jones, O. T., 1956. Quart. J. Geol. Soc. London, 111 (for 1955) 323. 

Jeffreys, H. & Bland, M. E. M., 1951. Mon. Not. R. Astr. Soc. Geophys. 
Suppl., 6, 148. 


OW 
ve | 
10- 
ve 
of | 
ng 
he | 
be 
on | 
ng | 
he | 
ly | 
er | 
m | 
ct, 
eS 
is : 
ht 
of 
he | 
is 
of | 
he | 
n. 
in | 
an 
on 
O- 
ne 
h 
th 
a 
y- 
an 
od 
es 
he | 
2d 
on 
at 
| 


360 Book Reviews 


Geophysics Series of the Bulletin of the Academy of Sciences of the 
U.S.S.R., No. 1. 


Pergamon Institute ({9. 10. 0. per year, monthly) 


This is a well-produced translation of what is presumably the Russian equiva- 
lent of the Geophysical Journal. Since geophysics is becoming more international 
in flavour and since many opportunities now occur for English-speaking and 
Russian geophysicists to meet, it is desirable that each should be able to read what 
progress the other is making. The Royal Astronomical Society has made a notable 
advance in one direction by providing summaries in Russian of all its published 
papers, and it is pleasing to see that the Pergamon Institute are now providing the 
full text of the Russian geophysical papers. 

The Russian journal makes a distinction between Terrestrial Physics and 
Atmospheric Physics. In the present number there are eight papers in the first 
group and three in the second. In this country the papers in the Atmospheric 
Physics group, which are about waterspouts, clouds and raindrops, would pro- 
bably have been in journals other than the Geophysical Journal, but all the papers 
on Terrestrial Physics are of the type that might be published in the Geophysical 
Journal. 

The first paper is a review article which describes the present position of seis- 
mology in the Chinese Peoples Republic and ends with a short fanfare praising 
the work of the Communist party. 

There are two seismological papers which relate observations to crustal struc- 
ture and a gravity paper giving a new method for reduction of results and cor- 
recting for topographical features. This last paper would really be more suitable 
for a publication dealing particularly with exploration geophysics, as would that 
on “Frequency Investigations in Detailed Electrical Prospecting”’. 

The wide variety of subjects that are covered includes magnetism, and results 
are given ot experimental studies of the stability of the magnetic fields of rocks. 
Another paper describes the equipment used to measure the natural electromag- 
netic field of the Earth in the frequency range 2-300c/s. This work has been done 
in order to develop a survey method. 

The last paper described the effect of turbulence on the formation of river- 
beds. It is concluded that the water currents of the Earth spend their energy largely 
in overcoming hydraulic resistance. and only to a few per cent on processes of 
erosion. 

There is a section of the journal which includes Reports and Notes on such 
practical subjects as gamma-ray well-logging and new seismograph design. T'wo 
of the Notes describe astronomical work (velocity of corpuscles projected from the 
Sun and the effect of meteoric streams on the electric field of the Earth’s atmos- 
phere) which would be placed by the R.A.S. in its Monthly Notices. 

The U.S.S.R. Geophysics Series No. 1 indicates that Russian thought is along 
the same lines as in this country. There is nothing particularly new in any of the 
papers in this first translation, and in some ways the journal is more suited to the 
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practical oil or mineral prospector than to the seeker after fundamental principles 
of the construction of the earth. ’ 

Some of the translation is a bit heavy; “expeditionary investigations” on p. 9 
and “‘cavernometry” on p. 130 sound a little odd. 

The price of the journal is really out of reach of the already overburdened 
individual geophysicist; he will be better served by buying the individual papers 
of interest to him from one of the several organizations that now provide a Russian 
translating service. The subjects covered are too diverse for any one research 
person but the journal might be useful to have in a reference library. 


T.F.G. 
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PESIOMEJJORJIAJIOB, B ITIEPEBOJIE HA PYCCKUM A3bIK 


CMEINMBAHME M OTJINBOB B Y83KUX MECTAX YCTbA 
PERM MEPCEM 


IT. 


Msamepenna B y3KuX MecTax peku Mepceii yKa3biBaloT Ha TO, YTO 
CMeIIMBaHHe H3-3a TYPOYeHTHOCTH BLISBAHHO! MO7KET OKOHYMTBCA Cepe MHI 
HbI€ CKOPOCTH, HallpaBseHHble B MOpe B BepXHeM C10e BBeEpX MO TeYeHMIO y ee 
3bIBaeT, YTO CHIbHO, TO NO MOE] He 
vty Teopulo K YCTbIO, B KOTOPOM 
peku (26 M3/cek.), BemunHy B 160 M2/cek. Bo Bpema 6omee 
CHJIbHOrO (103 M%/cek), Korga yc1oBuA MeHee 


O XUMHYECRKOM COCTABE HAPY#KHOPO AJIPA 3EMJIM 
Topdon Jlae. MarJ[onaavd u Jleonud Kuonoe 


CuntaetcA, 4TO CpeqHAA CTpykTypa 3emam copnagaeT CTpyKTypoii XOHAPUTHEIX Me- 
TeOpuToB. COBMeCTHMO C HacTOsel CKOPOCTbIO Terma 
B Bemae. CocTOAIee U3 HUKeIA, HEBO3MOAHO C XOH- 
MOeIbIO, KpoMe Toro oTHOWeHHe (Fe+Mg)/Si B o6 
COCTOAA M3 MaTepHada, lO CTpykType MMpoKceHy WIM 
B 400 KM, CTpyKTYpy KOTOPHIM CTPyKTypHI 
TO 10 MOfeTe ALPO COCTOATE U3 BeulecTBa (Fe, Ni),,¢Si. 

Brejen cpequuit Homep (Z) XapakTepHCTHKM yMepeHHo 
BLICOKUX Bprancnenue Z Ha OCHOBE OTHOLIe- 
HapymHoro spa Z nosyyaetca 22. nomep 
Fe,,.Si=22,8. Takum o6pasoM, BellleCTBO BbIBeeHHOe M3 ypaBHeHMA COCTOAHMA, 
COBMeCTMMO C Semun. 

IIpu CTpyKType Apa, (Fe, Ni);,.Si 
MecTHMa C TpeOoBaHHAMM TeOPHH MarHuTHOTO Semin. 


HOBOE PASHOCTH CWJIbl TATOTEHHA MEH] Y 
M3MEPEHWAMM B HAITMOHAJIBHOM JIABOPATOPHU B 
TIT IMHITTOHE B OCHOBHOM, A TAK}HKE BO BTOPMUHbIX CTAHIIMAX 

B AMPUKE 


Tog 


TuMa Ha OCHOBHOM OpuTaHCKOi CTaHIMM MO M3MepeHMIO TATOTeCHUA IPM 
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usMepeHuio TAroTeHuA B MoxannecOypre* Ha BTOPHYHEIX 
B IIperopuu u B Moy6pati. Bo spema MaxTHMKa HeywaAHHO 
noyBepriucb jelicTBuio pryTH. a uMeHHO IA, IB, IC VIA, onpeje- 
JIMIOT PasHOCTH CHJILI TTAOTeHMA 
= 1 ,5495+.0,0001 cm/cex* 
=2 ,6472+.0,0003 cm/cex? 
=2,5672+0,0003 cm/cek? 
B g=981,1963 cm/ceK? 
B Moy6pait g=979,6468 cm/ceK? 
B Moxannec6ypre g=978,5491 cm/cek? 
B IIperopun g=978,6291 cm/cek? 


Russian translations of summaries of original papers 


OTH BeIMYMHE C B 1948-9 rr. Xeiinca u Toda. 

*Iloka He pemieHo um cranyua B Moxannec6ypre cranunelt 
8B lOmnot Adpuke, HO OHA ABNACTCH 
ceTH. 


MATHUTHAH OUMCTKA TOPHbIX TOPO], TPM 


A. Acu JI. A. 


KOTOpOrO OOBIYHO MeCTHLIM SeMJIM B BpeMA 
BO-BTOPHIX, C CHJIbHOM CuO, KOTOpOrO 
BO BpeMA OOpa3z0BaHHA 

OnpeyeseHo, 4TO MO*KHO YCTPaHMTb 
HOM CHJILI, COXpaHMB 9TOM KOTMYeCTBO 
YCTHYHOTO OTO OCTMraeTCA CO“eTAHMeM MepeMeHHOTO BO3- 
M HarpeBaHHA. 

OGbIKHOBeHHO, C NepeMeHHBIMM 300 apcTey 
Temnepatype B 150°L]. [IA OUMCTKM TOPHEIX 


Y BEPETLA C.I.A. 32° 36° 
M 121° 128° 3.71. 


P. I. Meiicon 


IIpusegena Kapta MarHuTHOM MHTeHCHBHOCTH OOfacTH 32° u 36° 
121° 128° 3.4. Ha nett PA y3kux aHoMaJMii, 
B 400y , C CeBepa Ha lor Ha 500 KM CpefM 
AHOMAaJIMM Ha OKeaHa M COBMECTUMBI C 
4YMeM OTHOCHTebHO TOHKMX CJIO€B BelllecTBAa, 
MeHee Morin Ob BOBHUKHYTb M3-3a 
pesbeda BepxHelt MOBePXHOCTM MarHMTHOTO NOJIA, B 
C10EM OCAKOB. 

AnoMaJIMM MepeceKaloTcA 30HOM Myppoxv, ux OuepTaHMA Ta- 
06pa30M, 4YTO MO7KHO SoKOBOe BIpaBO Ha 
155 KM. 


KAJIMBPHPOBAHMA TIPMBOPOB JIJIH M3MEPEHMA 
TAVKECTU 


P. Tuavs6bepm 


u3MepeHHAX CHILI Ha 30HY M3MeHeHHA 
ee, BEIACHMJIOCh, NpHOOphl MeHee 4eM 
omm6oK, HalieHHEIX IPM TpeMA M3MepeHMA TA*KECTH, 
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MeCHAA a Tak*Ke pasOupaloTCA HEKOTOPLie BO3MO7KHHIe 
omM6oK. 4YTO MOCTOAHHEIA co BpemMeHeM Cy- 
IN€CTBEHHO M3MeHAWTCA M ONMCHIBACTCA CHCTeMa STO NOCTOAHHOM, BKIIOVeHHAA 
B AIA KOHTPOA NOCTOAHHOM B BpemaA. JlaHHEIe, 
npHOop OKasasICA B HeKOTOPHIX HO He HACTONbKO, 
HaCKOIbKO 

B BEIBOAMTCA, YTO TOUHOTO M3sMepeHHA CHIE 
TAKECTH, HEOOXOAMMO OCTOPO7KHO peryAAPHO MpoBepAA 
npuOopa u 4TO U3MepeHUA MOMOLIM MAATHMKOB B 
MaMepeHHit pasHOcTei TATOTeHMA. 


HABJIOTEHWUAM BTOPbIM CITYTHUKOM 3EMJIM (19578). 


A. X. Kyr. 


KOponeBCKOrO OOmjecTBAa CKOPOCTbIO BpaljeHuA 
BToporo CnyTHuKa coyeTawTcA ypaBHeHuAMM (1943 
1952) qua J u D, B SemiM 
no chepuyeckuM PyHKUMAM, C MepeoweHku J u D. HeO- 
pesyIbTATOB yKasbiBaeT HA WeHHOCTh Ha 


M Ha HeOOXOAMMOCTH CHILI TAKeECTH Ha 
Semsn. 


i 
ey 
4 
q 
ay 
a 
= 
= 
5 
. 


Geophysical Journal 


RECOMMENDATIONS TO AUTHORS 


1. Authors are advised to follow the general suggestions for the preparation of 
scientific papers given in Notes on the preparation of papers to be communicated to 
the Royal Society and in The printing of mathematics by T. W. Chaundy, P. R. 
Barrett and C. Batey (Oxford University Press, 1954). In particular the abbrevia- 
tions for the names of periodicals used in the Geophysical Fournal are those given 
in Notes on the preparation of papers . . . and the symbols, signs and abbreviations 


are those recommended for British scientific publications by the Symbols Commit- 
tee of the Royal Society. 


2. Papers should be typewritten with double spacing and on one side of the 
paper only. Two copies, one of them the top copy, should be submitted. All papers 


must be accompanied by a summary; summaries are not required for Letters to 
the Editors. 


3. Tables are printed without rules. They should be numbered serially with 
Arabic numerals. Table headings should be brief. Units should be placed at the 


head of the column. Tables should be typed on separate sheets and their positions 
in the text indicated on the copy. 


4. Illustrations, especially photographs, should be kept to a minimum, and the 
same information should not in general be given in both tables and illustrations. 
Line drawings should be in dense black ink on smooth white board or transparent 
tracing film. Original drawings must not be lettered; lettering should be indicated 
on copies or photoprints which must be-supplied in addition to the originals. 
Wherever possible, lettering should be kept outside the diagram so that it may be 
set in type. Diagrams should be drawn at twice to three times the size at which 
they will be printed. The maximum dimensions for diagrams, when reduced and 
with their lettering and captions, are 8 in by 5 in. A wide variety of flat tints of 
dots, lines and shadings can be applied to line drawings by the blockmaker. Photo- 


graphs for reproduction should be unmounted glossy prints and should be accom- 
panied by lettered prints. 


5. References should be quoted in the form recommended by the Royal Society 


but authors may also quote titles of papers if they wish. References in the text are 
made in the form (Smith 1940). 


6. Authors whose work involves much mathematics should follow carefully the 
recommendations in The Printing of Mathematics. 


7. The Editors will be glad to advise authors on any special points of difficulty 
arising in the preparation of papers for the Fournal. 
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